mini_project
Leung Chun Kit(21018713D), Lam Wai Tsun(21019252D)
Data source: https://archive.ics.uci.edu/ml/datasets/Beijing+Multi-Site+Air-Quality+Data
The data set includes hourly air pollutants data from 12 nationally-controlled air-quality monitoring sites. The air-quality data are from the Beijing Municipal Environmental Monitoring Center. The meteorological data in each air-quality site are matched with the nearest weather station from the China Meteorological Administration. The time period is from March 1st, 2013 to February 28th, 2017. Missing data are denoted as NA.
Attribute Information:
No: row number
year: year of data in this row
month: month of data in this row
day: day of data in this row
hour: hour of data in this row
PM2.5: PM2.5 concentration (ug/m^3)
PM10: PM10 concentration (ug/m^3)
SO2: SO2 concentration (ug/m^3)
NO2: NO2 concentration (ug/m^3)
CO: CO concentration (ug/m^3)
O3: O3 concentration (ug/m^3)
TEMP: temperature (degree Celsius)
PRES: pressure (hPa)
DEWP: dew point temperature (degree Celsius)
RAIN: precipitation (mm)
wd: wind direction
WSPM: wind speed (m/s)
station: name of the air-quality monitoring site
loading library
library(corrplot)
library(RColorBrewer)
library(tidyr)
library(dplyr)
library(data.table)
library(seasonal)
library(forecast)
library(dplyr)
library(TSstudio)
library(xts)
library(tseries)
library(ggplot2)
library(lattice)data <- read.csv("PRSA_Data_Wanshouxigong_20130301-20170228.csv", header = TRUE)We are going to study the air pollution condition in Beijing, Wanshou xigong.
attach(data)Check the number of missing value
sum(is.na(data))## [1] 5146fill the missing value
column_nane_list <- colnames(data)
for (a in column_nane_list){
data <-data %>% fill(a,.direction = 'updown')
}## Warning: Using an external vector in selections was deprecated in tidyselect 1.1.0.
## ℹ Please use `all_of()` or `any_of()` instead.
## # Was:
## data %>% select(a)
##
## # Now:
## data %>% select(all_of(a))
##
## See <https://tidyselect.r-lib.org/reference/faq-external-vector.html>.sum( apply(data,2, function(x) is.na(x)))## [1] 0We fill the missing value by the previous value.If the previous value is not available, we will fill by the next value.
There are two type of information in the dataset. One is the different index of the air pollutants (e.g. PM2.5, PM10, SO2, NO2, CO, O3). And the other type is the weather condition (e.g. TEMP, PRES, DEWP, RAIN, wd, WSPM).
grouped_year = aggregate(cbind(PM2.5,PM10,SO2,NO2,CO,O3,TEMP,PRES,DEWP,RAIN,WSPM) ~ data$year, data = data, FUN = mean, na.rm = TRUE)
grouped_month = aggregate(cbind(PM2.5,PM10,SO2,NO2,CO,O3,TEMP,PRES, DEWP,RAIN,WSPM) ~ data$year+ data$month, data = data, FUN = mean, na.rm = TRUE)
grouped_day = aggregate(cbind(PM2.5,PM10,SO2,NO2,CO,O3,TEMP,PRES, DEWP,RAIN,WSPM) ~ data$year+ data$month + data$day, data = data, FUN = mean, na.rm = TRUE)
grouped_year <- grouped_year[
order( grouped_year[,1] ),
]
grouped_month <- grouped_month[
order( grouped_month[,1], grouped_month[,2] ),
]
grouped_day <- grouped_day[
order( grouped_day[,1], grouped_day[,2], grouped_day[,3] ),
]Plot the linear correlation between diffrent variable
corr_month <- cor(grouped_month[,-c(1:2)])
corrplot(corr_month, type="upper", order="hclust",
col=brewer.pal(n=8, name="RdYlBu"))
We have found that the air pollutants have high correlation between each other (e.g. PM2.5, PM10, CO, SO2 and NO2 has positive correlation, but O3 has the negativity correlation with other air pollutants.)
And we also have found that the weather conditions also have high correlation with the air pollutants (e.g temperature, precipitation has negativity correlation with air pollutant(PM2.5,PM10, CO, SO2, NO2))
convert dataframe to timeframe object
grouped_month <- grouped_month [,-c(1:2)]
ts_by_month <- ts(grouped_month , start = c(2013, 3), frequency = 12)variables against time plot
plot.ts(ts_by_month[,1:6])
plot.ts(ts_by_month[,7:11])
# dates_by_month <- seq(as.Date("2013-03-01"),length = nrow(grouped_month), by ="month")
#
#
#
# xts_by_month <- xts(grouped_month, order.by = dates_by_month , frequency = 12)
#
# plot.xts(xts_by_month[,1:5], multi.panel = TRUE)
#
# plot.xts(xts_by_month[,11])
# ts_by_day <- xts(grouped_day, order.by = dates, frequency = 12)plot the box plot for variable in the cycle (month)
# cycle(ts_by_month)
boxplot(ts_by_month [,'PM2.5'] ~ cycle(ts_by_month [,'PM2.5']))
boxplot(ts_by_month [,'PM10'] ~ cycle(ts_by_month [,'PM10']))
boxplot(ts_by_month [,'SO2'] ~ cycle(ts_by_month [,'SO2']))
boxplot(ts_by_month [,'NO2'] ~ cycle(ts_by_month [,'NO2']))
boxplot(ts_by_month [,'CO'] ~ cycle(ts_by_month [,'CO']))
boxplot(ts_by_month [,'O3'] ~ cycle(ts_by_month [,'O3']))
boxplot(ts_by_month [,'TEMP'] ~ cycle(ts_by_month [,'TEMP']))
boxplot(ts_by_month [,'DEWP'] ~ cycle(ts_by_month [,'DEWP']))
boxplot(ts_by_month [,'PRES'] ~ cycle(ts_by_month [,'PRES']))
boxplot(ts_by_month [,'RAIN'] ~ cycle(ts_by_month [,'RAIN']))
boxplot(ts_by_month [,'WSPM'] ~ cycle(ts_by_month [,'WSPM']))
We can observe the distribution of difference variable among month
Plot the heat map
ts_heatmap(ts_by_month[,'PM2.5'],title = " Heatmap - the PM2.5 concentration in Wanshou xigong")ts_heatmap(ts_by_month[,'PM10'],title = " Heatmap - the PM10 concentration in Wanshou xigong")ts_heatmap(ts_by_month[,'SO2'],title = " Heatmap - the SO2 concentration in Wanshou xigong")ts_heatmap(ts_by_month[,'NO2'],title = " Heatmap - the NO2 concentration in Wanshou xigong")ts_heatmap(ts_by_month[,'CO'],title = " Heatmap - the CO concentration in Wanshou xigong")ts_heatmap(ts_by_month[,'O3'],title = " Heatmap - the O3 concentration in Wanshou xigong")ts_heatmap(ts_by_month[,'TEMP'],title = " Heatmap - the TEMP concentration in Wanshou xigong")ts_heatmap(ts_by_month[,'PRES'],title = " Heatmap - the PRES concentration in Wanshou xigong")ts_heatmap(ts_by_month[,'RAIN'],title = " Heatmap - the RAIN concentration in Wanshou xigong")Weekly Analysis
convert_date_to_weekday <- function(day, month, year){
date <- as.POSIXlt(paste(as.character(year) , as.character(month), as.character(day), sep="-"), tz = "UTC")
weekday <- weekdays(date)
return (weekday)
}data_weekday <- data
data_weekday$weekday <- mapply(convert_date_to_weekday, data$day, data$month, data$year)
grouped_weekday_year = aggregate(cbind(PM2.5,PM10,SO2,NO2,CO,O3,TEMP,PRES, DEWP,RAIN,WSPM) ~ year + weekday , data = data_weekday, FUN = mean, na.rm = TRUE)
grouped_weekday_year$weekday <- factor(grouped_weekday_year$weekday, levels = c("Monday", "Tuesday","Wednesday","Thursday", "Friday","Saturday","Sunday"))
weekdayorder = c("Monday", "Tuesday","Wednesday","Thursday", "Friday","Saturday","Sunday")
grouped_weekday_year = grouped_weekday_year[order(match(grouped_weekday_year$weekday,weekdayorder)),]
xyplot( PM2.5 ~ weekday ,
type="b",
group=year,
data = grouped_weekday_year,
auto.key = TRUE,
)
xyplot( CO ~ weekday ,
type="b",
group=year,
data = grouped_weekday_year,
auto.key = TRUE,
)
xyplot( TEMP ~ weekday ,
type="b",
group=year,
data = grouped_weekday_year,
auto.key = TRUE,
)
xyplot( NO2 ~ weekday ,
type="b",
group=year,
data = grouped_weekday_year,
auto.key = TRUE,
)
Hourly analysis
grouped_hour_year = aggregate(cbind(PM2.5,PM10,SO2,NO2,CO,O3,TEMP,PRES, DEWP,RAIN,WSPM) ~ year + hour , data = data, FUN = mean, na.rm = TRUE)
grouped_hour = aggregate(cbind(PM2.5,PM10,SO2,NO2,CO,O3,TEMP,PRES, DEWP,RAIN,WSPM) ~ hour , data = data, FUN = mean, na.rm = TRUE)
grouped_hour_median = aggregate(cbind(PM2.5,PM10,SO2,NO2,CO,O3,TEMP,PRES, DEWP,RAIN,WSPM) ~ hour , data = data, FUN = median, na.rm = TRUE)
library(lattice)
xyplot( PM2.5 ~ hour ,
type="b",
group=(year),
data = grouped_hour_year,
auto.key = TRUE,
)
xyplot( PM2.5 ~ hour ,
type="b",
data = grouped_hour,
auto.key = TRUE,
)
xyplot( PM2.5 ~ hour ,
type="b",
data = grouped_hour_median,
auto.key = TRUE,
)
xyplot( TEMP ~ hour ,
type="b",
group=(year),
data = grouped_hour_year,
auto.key = TRUE,
)
fit decomposition models (monthly) for PM2.5
# additive decomposition
add_decomp_PM25 <- decompose(ts_by_month[,'PM2.5'], type = 'additive')
add_decomp<- decompose(ts_by_month, type = 'additive')
add_decomp_PM25## $x
## Jan Feb Mar Apr May Jun Jul
## 2013 106.69892 78.82778 81.71640 105.67153 68.14651
## 2014 113.49194 156.06399 98.66667 86.64444 59.18804 58.03889 85.52151
## 2015 103.28199 102.95833 87.32258 71.85556 55.34543 59.44583 61.57124
## 2016 75.26210 48.35201 95.18952 67.05417 55.25806 64.08889 74.10753
## 2017 132.08333 77.61905
## Aug Sep Oct Nov Dec
## 2013 59.41532 61.74306 98.97849 84.30972 92.30914
## 2014 64.41331 70.99903 114.25806 103.31708 67.94247
## 2015 45.68548 53.15000 77.78226 125.42639 165.81989
## 2016 50.55242 59.25139 83.71640 106.37361 156.83468
## 2017
##
## $seasonal
## Jan Feb Mar Apr May Jun
## 2013 12.596727 -5.698226 -24.380240 -21.655511
## 2014 16.017806 21.170881 12.596727 -5.698226 -24.380240 -21.655511
## 2015 16.017806 21.170881 12.596727 -5.698226 -24.380240 -21.655511
## 2016 16.017806 21.170881 12.596727 -5.698226 -24.380240 -21.655511
## 2017 16.017806 21.170881
## Jul Aug Sep Oct Nov Dec
## 2013 -9.601027 -28.952748 -21.657639 13.707981 21.583772 26.868223
## 2014 -9.601027 -28.952748 -21.657639 13.707981 21.583772 26.868223
## 2015 -9.601027 -28.952748 -21.657639 13.707981 21.583772 26.868223
## 2016 -9.601027 -28.952748 -21.657639 13.707981 21.583772 26.868223
## 2017
##
## $trend
## Jan Feb Mar Apr May Jun Jul Aug
## 2013 NA NA NA NA NA NA
## 2014 87.14031 88.07252 88.66643 89.68875 91.11737 90.89406 89.45337 86.81522
## 2015 81.22384 79.44558 77.92155 75.65801 75.05941 80.05886 82.96959 79.52683
## 2016 78.40901 79.13415 79.59116 80.09264 79.54603 78.37778 80.37095 83.95796
## 2017 NA NA
## Sep Oct Nov Dec
## 2013 91.94639 91.93741 91.32442 88.40104
## 2014 84.12981 83.04094 82.26463 82.16314
## 2015 77.57935 77.70708 77.50339 77.69321
## 2016 NA NA NA NA
## 2017
##
## $random
## Jan Feb Mar Apr May
## 2013 NA NA NA
## 2014 10.33382068 46.82059124 -2.59649089 2.65392485 -7.54908995
## 2015 6.04034552 2.34186903 -3.19569409 1.89576779 4.66626068
## 2016 -19.16472204 -51.95301612 3.00162914 -7.34024848 0.09227343
## 2017 NA NA
## Jun Jul Aug Sep Oct
## 2013 NA NA NA -8.54569396 -6.66689181
## 2014 -11.19966299 5.66916168 6.55083331 8.52685202 17.50914281
## 2015 1.04248834 -11.79732274 -4.88859599 -2.77171390 -13.63280684
## 2016 7.36661881 3.33760522 -4.45279316 NA NA
## 2017
## Nov Dec
## 2013 -28.59846862 -22.96012658
## 2014 -0.53131766 -41.08889239
## 2015 26.33923044 61.25846313
## 2016 NA NA
## 2017
##
## $figure
## [1] 12.596727 -5.698226 -24.380240 -21.655511 -9.601027 -28.952748
## [7] -21.657639 13.707981 21.583772 26.868223 16.017806 21.170881
##
## $type
## [1] "additive"
##
## attr(,"class")
## [1] "decomposed.ts"plot(add_decomp_PM25)
add_decomp## $x
## PM2.5 PM10 SO2 NO2 CO O3
## Mar 2013 106.69892 130.19355 41.782387 63.42723 1563.6962 49.67568
## Apr 2013 78.82778 93.94167 22.014413 42.76404 1000.7417 50.78922
## May 2013 81.71640 129.03360 26.635918 52.47490 1094.1250 71.12242
## Jun 2013 105.67153 118.43681 14.815358 52.75863 1628.0139 69.63165
## Jul 2013 68.14651 73.75403 7.801311 45.31321 1223.7702 84.50378
## Aug 2013 59.41532 81.54301 6.996331 45.75652 1075.8065 86.94474
## Sep 2013 61.74306 117.49028 13.512820 57.73660 1326.5194 49.56147
## Oct 2013 98.97849 133.95833 15.550846 73.50200 1392.0699 28.22625
## Nov 2013 84.30972 111.67361 23.508333 58.00556 1575.4167 19.07793
## Dec 2013 92.30914 121.84812 44.816955 66.48069 2213.8441 12.87716
## Jan 2014 113.49194 152.32863 61.379032 79.61559 2171.3710 12.22043
## Feb 2014 156.06399 176.36310 67.059524 74.99405 2458.0357 18.06548
## Mar 2014 98.66667 157.95161 40.067204 66.39785 1489.1129 47.58333
## Apr 2014 86.64444 141.04306 18.758333 54.85417 1075.8333 68.79861
## May 2014 59.18804 112.96250 14.621640 47.10403 934.5430 92.11142
## Jun 2014 58.03889 81.15000 6.341667 46.40833 990.5556 103.63750
## Jul 2014 85.52151 110.45161 5.244624 40.52151 907.9301 104.85215
## Aug 2014 64.41331 91.18616 4.044624 49.18696 1027.1505 94.55497
## Sep 2014 70.99903 97.71708 5.679167 59.13194 1162.5000 53.78361
## Oct 2014 114.25806 143.46505 8.204301 76.72177 1372.3118 26.74059
## Nov 2014 103.31708 142.29014 14.245139 77.05292 1785.4167 22.32083
## Dec 2014 67.94247 111.75591 33.370027 63.90403 1873.9247 29.50847
## Jan 2015 103.28199 128.48925 39.075672 76.75376 2011.6935 24.22755
## Feb 2015 102.95833 133.78869 29.120536 59.45833 1557.1429 38.70536
## Mar 2015 87.32258 155.15323 23.263441 59.77419 1257.3925 54.28763
## Apr 2015 71.85556 121.17083 9.693056 49.62222 835.8333 78.03472
## May 2015 55.34543 93.78898 9.059140 41.90726 779.1667 98.06855
## Jun 2015 59.44583 81.00139 6.877778 37.44861 865.8333 92.75417
## Jul 2015 61.57124 77.40188 4.935484 38.15995 978.3602 95.05780
## Aug 2015 45.68548 67.98118 3.108871 35.40726 909.6774 100.36694
## Sep 2015 53.15000 68.63472 4.965278 43.25139 956.6667 58.35278
## Oct 2015 77.78226 98.78629 5.424731 55.75672 1043.9516 36.92473
## Nov 2015 125.42639 131.26806 14.548611 60.19444 2085.4167 14.15139
## Dec 2015 165.81989 175.36425 24.176075 73.86290 3024.3280 17.96102
## Jan 2016 75.26210 88.44892 24.372312 52.94892 1639.1129 28.16935
## Feb 2016 48.35201 61.49856 17.149425 34.41810 933.9080 47.92385
## Mar 2016 95.18952 130.83737 20.548387 53.99194 1223.9247 48.20565
## Apr 2016 67.05417 109.42500 11.688889 41.57083 852.0833 72.54167
## May 2016 55.25806 95.65995 9.151882 37.45833 724.3280 99.22043
## Jun 2016 64.08889 82.17778 5.462500 34.99444 865.8333 126.60556
## Jul 2016 74.10753 82.06452 3.322581 34.48925 997.9839 96.86290
## Aug 2016 50.55242 64.67608 2.659946 39.72849 932.6613 76.13575
## Sep 2016 59.25139 70.65417 2.827778 50.75833 870.1389 56.39444
## Oct 2016 83.71640 94.86962 3.310484 61.01210 1246.5054 20.25538
## Nov 2016 106.37361 136.37778 10.468056 71.47778 1812.7778 14.05556
## Dec 2016 156.83468 169.18011 17.418011 89.91532 2604.1667 15.15188
## Jan 2017 132.08333 149.47984 18.594086 71.23522 2254.3011 27.93952
## Feb 2017 77.61905 90.79464 18.520833 56.49405 1173.6607 40.98810
## TEMP PRES DEWP RAIN WSPM
## Mar 2013 6.373252688 1011.6395 -6.3188172 0.0236559140 1.753495
## Apr 2013 12.599166667 1007.2869 -2.8651389 0.0141666667 2.223472
## May 2013 21.787903226 1002.1239 8.4520161 0.0076612903 1.771237
## Jun 2013 23.783055556 999.0772 17.5343056 0.1015277778 1.355000
## Jul 2013 27.156451613 995.0235 21.1831989 0.2802419355 1.359946
## Aug 2013 27.257661290 998.2694 20.0935484 0.0907258065 1.416398
## Sep 2013 20.487361111 1008.3138 14.5504167 0.1150000000 1.160972
## Oct 2013 13.199865591 1015.2367 5.3115591 0.0170698925 1.337366
## Nov 2013 5.880555556 1016.0646 -7.1959722 0.0008333333 1.769028
## Dec 2013 -0.005241935 1019.1534 -13.2763441 0.0000000000 1.699866
## Jan 2014 0.173655914 1019.2828 -12.6904570 0.0000000000 1.618548
## Feb 2014 0.089434524 1021.6240 -9.5802083 0.0098214286 1.438244
## Mar 2014 10.513844086 1012.8694 -5.8724462 0.0001344086 1.810215
## Apr 2014 17.184027778 1009.5782 3.9636111 0.0180555556 1.615556
## May 2014 21.865188172 1000.7228 6.9375000 0.0915322581 2.020027
## Jun 2014 24.901944444 999.3411 16.3420833 0.1980555556 1.393472
## Jul 2014 28.287634409 997.8254 19.5525538 0.0581989247 1.479973
## Aug 2014 26.208736559 1001.8782 18.1295699 0.1381720430 1.395833
## Sep 2014 21.123472222 1007.9142 13.8113889 0.0547222222 1.310833
## Oct 2014 14.044892473 1013.6585 6.2706989 0.0182795699 1.217070
## Nov 2014 6.371111111 1021.8025 -5.1972222 0.0004166667 1.538472
## Dec 2014 -0.226478495 1025.9004 -16.1875000 0.0000000000 2.189113
## Jan 2015 -0.125537634 1020.9488 -13.7850806 0.0005376344 1.604704
## Feb 2015 1.743154762 1017.6015 -12.4601190 0.0168154762 1.811905
## Mar 2015 8.439861751 1016.8734 -9.8809140 0.0088709677 2.213575
## Apr 2015 15.630138889 1010.8733 1.6945833 0.0556944444 2.483056
## May 2015 21.455510753 1003.4489 6.6134409 0.0501344086 2.471102
## Jun 2015 24.537777778 1000.3461 14.3013889 0.1362500000 1.980278
## Jul 2015 26.397849462 1001.2784 18.5057796 0.2701612903 1.627285
## Aug 2015 26.435752688 1003.4042 18.4461022 0.1258064516 1.615457
## Sep 2015 20.600833333 1011.3549 13.7698611 0.1372222222 1.601944
## Oct 2015 14.319758065 1015.7254 4.1530914 0.0190860215 1.848790
## Nov 2015 2.847777778 1023.4250 -1.2023611 0.0522222222 1.548056
## Dec 2015 -0.362768817 1024.0706 -7.2271505 0.0028225806 1.856452
## Jan 2016 -4.550000000 1025.9500 -18.0596774 0.0006720430 2.115054
## Feb 2016 1.398419540 1022.5598 -15.0882184 0.0181034483 2.449138
## Mar 2016 9.201881720 1016.7401 -8.7526882 0.0000000000 2.049866
## Apr 2016 16.454722222 1007.8117 -0.3329167 0.0066666667 2.382361
## May 2016 21.672849462 1005.5414 6.2018817 0.0458333333 2.207796
## Jun 2016 25.735833333 1000.7619 14.6834722 0.0991666667 1.803056
## Jul 2016 27.699462366 1000.6353 20.7272849 0.4233870968 1.698925
## Aug 2016 27.863306452 1003.5901 19.2868280 0.0724462366 1.618414
## Sep 2016 22.197344026 1009.3406 14.0559722 0.1598611111 1.554444
## Oct 2016 13.338648073 1017.1267 7.2780914 0.1174731183 1.476882
## Nov 2016 4.414027778 1020.5669 -3.7584722 0.0073611111 1.557778
## Dec 2016 0.640725806 1023.4661 -9.2724462 0.0000000000 1.435618
## Jan 2017 -1.140412186 1025.0380 -12.7756720 0.0004032258 1.969220
## Feb 2017 2.661532738 1022.0206 -13.6276786 0.0061011905 1.953869
##
## $seasonal
## Jan Feb Mar Apr May Jun
## 2013 3.998767 -34.759863 -46.821312 -39.493936
## 2014 55.312609 29.866751 3.998767 -34.759863 -46.821312 -39.493936
## 2015 55.312609 29.866751 3.998767 -34.759863 -46.821312 -39.493936
## 2016 55.312609 29.866751 3.998767 -34.759863 -46.821312 -39.493936
## 2017 55.312609 29.866751 3.998767 -34.759863 -46.821312 -39.493936
## 2018 55.312609 29.866751 3.998767 -34.759863 -46.821312 -39.493936
## 2019 55.312609 29.866751 3.998767 -34.759863 -46.821312 -39.493936
## 2020 55.312609 29.866751 3.998767 -34.759863 -46.821312 -39.493936
## 2021 55.312609 29.866751 3.998767 -34.759863 -46.821312 -39.493936
## 2022 55.312609 29.866751 3.998767 -34.759863 -46.821312 -39.493936
## 2023 55.312609 29.866751 3.998767 -34.759863 -46.821312 -39.493936
## 2024 55.312609 29.866751 3.998767 -34.759863 -46.821312 -39.493936
## 2025 55.312609 29.866751 3.998767 -34.759863 -46.821312 -39.493936
## 2026 55.312609 29.866751 3.998767 -34.759863 -46.821312 -39.493936
## 2027 55.312609 29.866751 3.998767 -34.759863 -46.821312 -39.493936
## 2028 55.312609 29.866751 3.998767 -34.759863 -46.821312 -39.493936
## 2029 55.312609 29.866751 3.998767 -34.759863 -46.821312 -39.493936
## 2030 55.312609 29.866751 3.998767 -34.759863 -46.821312 -39.493936
## 2031 55.312609 29.866751 3.998767 -34.759863 -46.821312 -39.493936
## 2032 55.312609 29.866751 3.998767 -34.759863 -46.821312 -39.493936
## 2033 55.312609 29.866751 3.998767 -34.759863 -46.821312 -39.493936
## 2034 55.312609 29.866751 3.998767 -34.759863 -46.821312 -39.493936
## 2035 55.312609 29.866751 3.998767 -34.759863 -46.821312 -39.493936
## 2036 55.312609 29.866751 3.998767 -34.759863 -46.821312 -39.493936
## 2037 55.312609 29.866751 3.998767 -34.759863 -46.821312 -39.493936
## 2038 55.312609 29.866751 3.998767 -34.759863 -46.821312 -39.493936
## 2039 55.312609 29.866751 3.998767 -34.759863 -46.821312 -39.493936
## 2040 55.312609 29.866751 3.998767 -34.759863 -46.821312 -39.493936
## 2041 55.312609 29.866751 3.998767 -34.759863 -46.821312 -39.493936
## 2042 55.312609 29.866751 3.998767 -34.759863 -46.821312 -39.493936
## 2043 55.312609 29.866751 3.998767 -34.759863 -46.821312 -39.493936
## 2044 55.312609 29.866751 3.998767 -34.759863 -46.821312 -39.493936
## 2045 55.312609 29.866751 3.998767 -34.759863 -46.821312 -39.493936
## 2046 55.312609 29.866751 3.998767 -34.759863 -46.821312 -39.493936
## 2047 55.312609 29.866751 3.998767 -34.759863 -46.821312 -39.493936
## 2048 55.312609 29.866751 3.998767 -34.759863 -46.821312 -39.493936
## 2049 55.312609 29.866751 3.998767 -34.759863 -46.821312 -39.493936
## 2050 55.312609 29.866751 3.998767 -34.759863 -46.821312 -39.493936
## 2051 55.312609 29.866751 3.998767 -34.759863 -46.821312 -39.493936
## 2052 55.312609 29.866751 3.998767 -34.759863 -46.821312 -39.493936
## 2053 55.312609 29.866751 3.998767 -34.759863 -46.821312 -39.493936
## 2054 55.312609 29.866751 3.998767 -34.759863 -46.821312 -39.493936
## 2055 55.312609 29.866751 3.998767 -34.759863 -46.821312 -39.493936
## 2056 55.312609 29.866751 3.998767 -34.759863 -46.821312 -39.493936
## 2057 55.312609 29.866751
## Jul Aug Sep Oct Nov Dec
## 2013 -33.856298 -36.175662 -24.854130 -8.723305 40.661334 94.845043
## 2014 -33.856298 -36.175662 -24.854130 -8.723305 40.661334 94.845043
## 2015 -33.856298 -36.175662 -24.854130 -8.723305 40.661334 94.845043
## 2016 -33.856298 -36.175662 -24.854130 -8.723305 40.661334 94.845043
## 2017 -33.856298 -36.175662 -24.854130 -8.723305 40.661334 94.845043
## 2018 -33.856298 -36.175662 -24.854130 -8.723305 40.661334 94.845043
## 2019 -33.856298 -36.175662 -24.854130 -8.723305 40.661334 94.845043
## 2020 -33.856298 -36.175662 -24.854130 -8.723305 40.661334 94.845043
## 2021 -33.856298 -36.175662 -24.854130 -8.723305 40.661334 94.845043
## 2022 -33.856298 -36.175662 -24.854130 -8.723305 40.661334 94.845043
## 2023 -33.856298 -36.175662 -24.854130 -8.723305 40.661334 94.845043
## 2024 -33.856298 -36.175662 -24.854130 -8.723305 40.661334 94.845043
## 2025 -33.856298 -36.175662 -24.854130 -8.723305 40.661334 94.845043
## 2026 -33.856298 -36.175662 -24.854130 -8.723305 40.661334 94.845043
## 2027 -33.856298 -36.175662 -24.854130 -8.723305 40.661334 94.845043
## 2028 -33.856298 -36.175662 -24.854130 -8.723305 40.661334 94.845043
## 2029 -33.856298 -36.175662 -24.854130 -8.723305 40.661334 94.845043
## 2030 -33.856298 -36.175662 -24.854130 -8.723305 40.661334 94.845043
## 2031 -33.856298 -36.175662 -24.854130 -8.723305 40.661334 94.845043
## 2032 -33.856298 -36.175662 -24.854130 -8.723305 40.661334 94.845043
## 2033 -33.856298 -36.175662 -24.854130 -8.723305 40.661334 94.845043
## 2034 -33.856298 -36.175662 -24.854130 -8.723305 40.661334 94.845043
## 2035 -33.856298 -36.175662 -24.854130 -8.723305 40.661334 94.845043
## 2036 -33.856298 -36.175662 -24.854130 -8.723305 40.661334 94.845043
## 2037 -33.856298 -36.175662 -24.854130 -8.723305 40.661334 94.845043
## 2038 -33.856298 -36.175662 -24.854130 -8.723305 40.661334 94.845043
## 2039 -33.856298 -36.175662 -24.854130 -8.723305 40.661334 94.845043
## 2040 -33.856298 -36.175662 -24.854130 -8.723305 40.661334 94.845043
## 2041 -33.856298 -36.175662 -24.854130 -8.723305 40.661334 94.845043
## 2042 -33.856298 -36.175662 -24.854130 -8.723305 40.661334 94.845043
## 2043 -33.856298 -36.175662 -24.854130 -8.723305 40.661334 94.845043
## 2044 -33.856298 -36.175662 -24.854130 -8.723305 40.661334 94.845043
## 2045 -33.856298 -36.175662 -24.854130 -8.723305 40.661334 94.845043
## 2046 -33.856298 -36.175662 -24.854130 -8.723305 40.661334 94.845043
## 2047 -33.856298 -36.175662 -24.854130 -8.723305 40.661334 94.845043
## 2048 -33.856298 -36.175662 -24.854130 -8.723305 40.661334 94.845043
## 2049 -33.856298 -36.175662 -24.854130 -8.723305 40.661334 94.845043
## 2050 -33.856298 -36.175662 -24.854130 -8.723305 40.661334 94.845043
## 2051 -33.856298 -36.175662 -24.854130 -8.723305 40.661334 94.845043
## 2052 -33.856298 -36.175662 -24.854130 -8.723305 40.661334 94.845043
## 2053 -33.856298 -36.175662 -24.854130 -8.723305 40.661334 94.845043
## 2054 -33.856298 -36.175662 -24.854130 -8.723305 40.661334 94.845043
## 2055 -33.856298 -36.175662 -24.854130 -8.723305 40.661334 94.845043
## 2056 -33.856298 -36.175662 -24.854130 -8.723305 40.661334 94.845043
## 2057
##
## $trend
## [,1] [,2] [,3] [,4] [,5] [,6] [,7]
## Mar 2013 NA NA NA NA NA NA NA
## Apr 2013 NA NA NA NA NA NA NA
## May 2013 NA NA NA NA NA NA NA
## Jun 2013 NA NA NA NA NA NA NA
## Jul 2013 NA NA NA NA NA NA NA
## Aug 2013 NA NA NA NA NA NA NA
## Sep 2013 91.94639 121.20365 28.75130 59.52619 1557.177 45.97084 13.40445
## Oct 2013 91.93741 124.32279 28.54417 60.15372 1557.198 46.63405 13.76801
## Nov 2013 91.32442 125.61572 27.90790 60.43369 1553.677 48.25898 13.96227
## Dec 2013 88.40104 123.39247 27.05424 59.94531 1520.467 50.55043 14.01211
## Jan 2014 87.14031 123.36792 26.59464 59.48106 1480.747 52.81519 14.10586
## Feb 2014 88.07252 125.29879 26.36512 59.42434 1465.559 53.98013 14.10929
## Mar 2014 88.66643 124.87670 25.91573 59.62542 1456.698 54.47315 14.09209
## Apr 2014 89.68875 124.44893 25.28322 59.81771 1449.040 54.58717 14.15380
## May 2014 91.11737 126.12073 24.59115 60.74551 1456.967 54.66039 14.20945
## Jun 2014 90.89406 126.97591 23.72823 61.43179 1451.554 55.48848 14.22067
## Jul 2014 89.45337 125.56210 22.32197 61.20519 1430.737 56.68175 14.19899
## Aug 2014 86.81522 122.79486 19.81187 60.43862 1386.547 58.04204 14.25543
## Sep 2014 84.12981 120.90432 17.53092 59.51532 1339.355 59.18138 14.23792
## Oct 2014 83.04094 119.95971 16.45304 59.02133 1319.700 59.84556 14.08676
## Nov 2014 82.26463 118.33281 15.84355 58.58680 1303.226 60.47862 14.00494
## Dec 2014 82.16314 117.52772 15.63412 57.99695 1291.555 60.27336 13.97270
## Jan 2015 81.22384 116.14445 15.64358 57.52523 1289.293 59.41179 13.87878
## Feb 2015 79.44558 113.80051 15.59171 56.85268 1287.332 59.24585 13.80950
## Mar 2015 77.92155 111.62187 15.52297 55.61683 1273.861 59.67840 13.79718
## Apr 2015 75.65801 108.54849 15.37741 54.08160 1251.603 60.29312 13.78686
## May 2015 75.05941 106.22762 15.27424 52.50562 1250.422 60.37707 13.65151
## Jun 2015 80.05886 108.41871 14.90380 52.21813 1310.855 59.55553 13.49902
## Jul 2015 82.96959 109.40072 13.90808 51.64122 1343.264 59.23863 13.30899
## Aug 2015 79.52683 104.72028 12.79665 49.60601 1301.772 59.78697 13.11027
## Sep 2015 77.57935 100.69503 12.18472 48.32174 1274.409 59.91766 13.12766
## Oct 2015 77.70708 99.19246 12.15476 47.74534 1273.692 59.43537 13.19377
## Nov 2015 77.50339 98.78101 12.24178 47.22449 1272.084 59.25449 13.23718
## Dec 2015 77.69321 98.90798 12.18667 46.93686 1269.799 60.71295 13.29616
## Jan 2016 78.40901 99.15127 12.06050 46.68166 1270.617 62.19864 13.40031
## Feb 2016 79.13415 99.20784 11.97459 46.70876 1272.392 61.26422 13.51402
## Mar 2016 79.59116 99.15427 11.86682 47.20160 1269.745 60.17299 13.64003
## Apr 2016 80.09264 99.07522 11.68967 47.73337 1274.579 59.39684 13.66567
## May 2016 79.54603 99.12493 11.43155 48.42248 1271.659 58.69829 13.69005
## Jun 2016 78.37778 99.08016 10.97994 49.56147 1242.792 58.57725 13.79712
## Jul 2016 80.37095 101.36544 10.45759 50.99225 1250.918 58.45062 13.98100
## Aug 2016 83.95796 105.12907 10.27398 52.67401 1286.541 58.15206 14.17570
## Sep 2016 NA NA NA NA NA NA NA
## Oct 2016 NA NA NA NA NA NA NA
## Nov 2016 NA NA NA NA NA NA NA
## Dec 2016 NA NA NA NA NA NA NA
## Jan 2017 NA NA NA NA NA NA NA
## Feb 2017 NA NA NA NA NA NA NA
## [,8] [,9] [,10] [,11]
## Mar 2013 NA NA NA NA
## Apr 2013 NA NA NA NA
## May 2013 NA NA NA NA
## Jun 2013 NA NA NA NA
## Jul 2013 NA NA NA NA
## Aug 2013 NA NA NA NA
## Sep 2013 1009.476 2.951774 0.05407861 1.577661
## Oct 2013 1009.623 3.254904 0.05326058 1.554694
## Nov 2013 1009.660 3.476331 0.05691724 1.539731
## Dec 2013 1009.612 3.363550 0.06443386 1.551700
## Jan 2014 1009.740 3.245931 0.05920406 1.558304
## Feb 2014 1010.007 3.096155 0.05192919 1.562449
## Mar 2014 1010.141 2.983529 0.05139454 1.567836
## Apr 2014 1010.058 2.992701 0.04893337 1.569068
## May 2014 1010.232 3.115946 0.04896641 1.554449
## Jun 2014 1010.752 3.077929 0.04894905 1.565228
## Jul 2014 1011.103 2.911022 0.04897145 1.585036
## Aug 2014 1011.004 2.745416 0.04928527 1.600029
## Sep 2014 1011.004 2.458400 0.04994072 1.632404
## Oct 2014 1011.224 2.196838 0.05187303 1.685357
## Nov 2014 1011.392 2.088793 0.05171640 1.740298
## Dec 2014 1011.547 1.990261 0.04741626 1.783543
## Jan 2015 1011.733 1.861617 0.05367280 1.814131
## Feb 2015 1011.941 1.831190 0.06198933 1.829420
## Mar 2015 1012.148 1.842649 0.06491159 1.850700
## Apr 2015 1012.377 1.752685 0.06838270 1.889152
## May 2015 1012.531 1.830903 0.07057486 1.915873
## Jun 2015 1012.522 2.370704 0.07285104 1.902411
## Jul 2015 1012.654 2.565944 0.07297424 1.909815
## Aug 2015 1013.069 2.278331 0.07303351 1.957631
## Sep 2015 1013.270 2.215837 0.07271755 1.977361
## Oct 2015 1013.137 2.178367 0.07030510 1.966344
## Nov 2015 1013.097 2.076739 0.06808307 1.951177
## Dec 2015 1013.201 2.075511 0.06635872 1.932822
## Jan 2016 1013.192 2.183994 0.07119799 1.928423
## Feb 2016 1013.173 2.311587 0.07535905 1.931531
## Mar 2016 1013.097 2.358539 0.07407900 1.929675
## Apr 2016 1013.071 2.500668 0.07912175 1.912200
## May 2016 1013.010 2.524372 0.08135200 1.897109
## Jun 2016 1012.866 2.332647 0.07936518 1.879979
## Jul 2016 1012.803 2.467593 0.07923637 1.856368
## Aug 2016 1012.742 2.748616 0.07872507 1.829655
## Sep 2016 NA NA NA NA
## Oct 2016 NA NA NA NA
## Nov 2016 NA NA NA NA
## Dec 2016 NA NA NA NA
## Jan 2017 NA NA NA NA
## Feb 2017 NA NA NA NA
##
## $random
## x - seasonal.x.PM2.5 x - seasonal.x.PM10 x - seasonal.x.SO2
## Mar 2013 NA NA NA
## Apr 2013 NA NA NA
## May 2013 NA NA NA
## Jun 2013 NA NA NA
## Jul 2013 NA NA NA
## Aug 2013 NA NA NA
## Sep 2013 -5.349203 21.1407608 9.6156468
## Oct 2013 15.764394 18.3588480 -4.2700155
## Nov 2013 -47.676030 -54.6034424 -45.0609031
## Dec 2013 -90.936947 -96.3893983 -77.0823253
## Jan 2014 -28.960982 -26.3519018 -20.5282141
## Feb 2014 38.124721 21.1975582 10.8276512
## Mar 2014 6.001469 29.0761455 10.1527061
## Apr 2014 31.715562 51.3539876 28.2349730
## May 2014 14.891982 33.6630789 36.8518008
## Jun 2014 6.638761 -6.3319773 22.1073734
## Jul 2014 29.924432 18.7458137 16.7789545
## Aug 2014 13.773747 4.5669619 20.4084164
## Sep 2014 11.723343 1.6668906 13.0023757
## Oct 2014 39.940429 32.2286453 0.4745623
## Nov 2014 -19.608879 -16.7040032 -42.2597487
## Dec 2014 -109.065713 -100.6168484 -77.1091370
## Jan 2015 -33.254457 -42.9678161 -31.8805143
## Feb 2015 -6.354002 -9.8785696 -16.3379228
## Mar 2015 5.402266 39.5325890 3.7417018
## Apr 2015 30.957405 47.3822063 29.0755069
## May 2015 27.107332 34.3826690 40.6062107
## Jun 2015 18.880913 12.0766096 31.4679093
## Jul 2015 12.457948 1.8574640 24.8836988
## Aug 2015 2.334318 -0.5634354 26.4878862
## Sep 2015 0.424777 -7.2061782 17.6346848
## Oct 2015 8.798479 8.3171360 1.9932813
## Nov 2015 7.261669 -8.1742857 -38.3545024
## Dec 2015 -6.718357 -18.3887762 -82.8556417
## Jan 2016 -58.459524 -66.0149568 -43.0007962
## Feb 2016 -60.648887 -67.5760246 -24.6919160
## Mar 2016 11.599589 27.6843314 4.6827979
## Apr 2016 21.721389 45.1096464 34.7590854
## May 2016 22.533345 43.3563308 44.5416439
## Jun 2016 25.205043 22.5915538 33.9764954
## Jul 2016 27.592876 14.5553719 26.7212835
## Aug 2016 2.770121 -4.2773296 28.5616305
## Sep 2016 NA NA NA
## Oct 2016 NA NA NA
## Nov 2016 NA NA NA
## Dec 2016 NA NA NA
## Jan 2017 NA NA NA
## Feb 2017 NA NA NA
## x - seasonal.x.NO2 x - seasonal.x.CO x - seasonal.x.O3
## Mar 2013 NA NA NA
## Apr 2013 NA NA NA
## May 2013 NA NA NA
## Jun 2013 NA NA NA
## Jul 2013 NA NA NA
## Aug 2013 NA NA NA
## Sep 2013 23.0645354 -205.80297 28.444759
## Oct 2013 22.0715805 -156.40453 -9.684491
## Nov 2013 -43.0894713 -18.92196 -69.842382
## Dec 2013 -88.3096672 598.53176 -132.518312
## Jan 2014 -35.1780781 635.31185 -95.907370
## Feb 2014 -14.2970456 962.60978 -65.781408
## Mar 2014 2.7736658 28.41643 -10.888582
## Apr 2014 29.7963159 -338.44712 48.971305
## May 2014 33.1798330 -475.60274 84.272349
## Jun 2014 24.4704786 -421.50426 87.642956
## Jul 2014 13.1726161 -488.95082 82.026701
## Aug 2014 24.9240014 -323.22059 72.688596
## Sep 2014 24.4707589 -152.00044 19.456361
## Oct 2014 26.4237475 61.33558 -24.381667
## Nov 2014 -22.1952197 441.52979 -78.819116
## Dec 2014 -88.9379592 487.52492 -125.609932
## Jan 2015 -36.0840733 667.08835 -90.496841
## Feb 2015 -27.2610934 239.94364 -50.407247
## Mar 2015 0.1585953 -20.46766 -9.389533
## Apr 2015 30.3004873 -381.01011 52.501464
## May 2015 36.2229523 -424.43365 84.512793
## Jun 2015 24.7244127 -405.52783 72.692573
## Jul 2015 20.3750252 -331.04786 69.675465
## Aug 2015 21.9769122 -355.91898 76.755624
## Sep 2015 19.7837817 -292.88866 23.289247
## Oct 2015 16.7346907 -221.01713 -13.787331
## Nov 2015 -27.6913785 772.67115 -85.764431
## Dec 2015 -67.9190000 1659.68368 -137.596977
## Jan 2016 -49.0453411 313.18341 -89.341896
## Feb 2016 -42.1574108 -368.35091 -43.207123
## Mar 2016 2.7915647 -49.81858 -15.966114
## Apr 2016 28.5973288 -387.73576 47.904690
## May 2016 35.8571648 -500.50947 87.343452
## Jun 2016 24.9269102 -337.46480 107.522242
## Jul 2016 17.3532955 -219.07802 72.268576
## Aug 2016 23.2301471 -317.70377 54.159356
## Sep 2016 NA NA NA
## Oct 2016 NA NA NA
## Nov 2016 NA NA NA
## Dec 2016 NA NA NA
## Jan 2017 NA NA NA
## Feb 2017 NA NA NA
## x - seasonal.x.TEMP x - seasonal.x.PRES x - seasonal.x.DEWP
## Mar 2013 NA NA NA
## Apr 2013 NA NA NA
## May 2013 NA NA NA
## Jun 2013 NA NA NA
## Jul 2013 NA NA NA
## Aug 2013 NA NA NA
## Sep 2013 31.937040 23.6920012 36.45277
## Oct 2013 8.155159 14.3374083 10.77996
## Nov 2013 -48.743047 -34.2564322 -51.33364
## Dec 2013 -108.862394 -85.3039816 -111.48494
## Jan 2014 -69.244814 -45.7698519 -71.24900
## Feb 2014 -43.886606 -18.2499467 -42.54311
## Mar 2014 -7.577011 -1.2702862 -12.85474
## Apr 2014 37.790088 34.2795918 35.73077
## May 2014 54.477048 37.3123752 50.64287
## Jun 2014 50.175207 28.0830542 52.75809
## Jul 2014 47.944943 20.5791651 50.49783
## Aug 2014 48.128971 27.0495383 51.55982
## Sep 2014 31.739686 21.7647147 36.20712
## Oct 2014 8.681443 11.1573916 12.79717
## Nov 2014 -48.295163 -30.2507663 -47.94735
## Dec 2014 -109.044218 -80.4920337 -113.02280
## Jan 2015 -69.316928 -46.0969601 -70.95931
## Feb 2015 -41.933096 -24.2058591 -44.15806
## Mar 2015 -9.356087 0.7271035 -15.72233
## Apr 2015 36.603143 33.2561724 34.70176
## May 2015 54.625317 37.7394862 51.60385
## Jun 2015 50.532692 27.3179356 51.42462
## Jul 2015 46.945157 22.4804063 49.79613
## Aug 2015 49.501141 26.5105980 52.34343
## Sep 2015 32.327303 22.9387221 36.40815
## Oct 2015 9.849295 11.3115657 10.69803
## Nov 2015 -51.050738 -30.3330940 -43.94043
## Dec 2015 -108.503969 -83.9757514 -104.14771
## Jan 2016 -73.262918 -42.5544155 -75.55628
## Feb 2016 -41.982357 -20.4797413 -47.26656
## Mar 2016 -8.436913 -0.3552892 -15.10999
## Apr 2016 37.548916 29.5004957 31.92628
## May 2016 54.804111 39.3523744 50.49882
## Jun 2016 51.432646 27.3898155 51.84476
## Jul 2016 47.574759 21.6887686 52.11599
## Aug 2016 49.863271 27.0233042 52.71387
## Sep 2016 NA NA NA
## Oct 2016 NA NA NA
## Nov 2016 NA NA NA
## Dec 2016 NA NA NA
## Jan 2017 NA NA NA
## Feb 2017 NA NA NA
## x - seasonal.x.RAIN x - seasonal.x.WSPM
## Mar 2013 NA NA
## Apr 2013 NA NA
## May 2013 NA NA
## Jun 2013 NA NA
## Jul 2013 NA NA
## Aug 2013 NA NA
## Sep 2013 24.915051 24.437441
## Oct 2013 8.687115 8.505977
## Nov 2013 -40.717418 -40.432037
## Dec 2013 -94.909477 -94.696878
## Jan 2014 -55.371813 -55.252364
## Feb 2014 -29.908859 -29.990956
## Mar 2014 -4.050027 -3.756388
## Apr 2014 34.728985 34.806350
## May 2014 46.863877 47.286889
## Jun 2014 39.643042 39.322180
## Jul 2014 33.865525 33.751234
## Aug 2014 36.264548 35.971466
## Sep 2014 24.858911 24.532559
## Oct 2014 8.689712 8.255018
## Nov 2014 -40.712634 -40.863160
## Dec 2014 -94.892460 -94.439473
## Jan 2015 -55.365744 -55.522035
## Feb 2015 -29.911925 -29.884266
## Mar 2015 -4.054808 -3.635892
## Apr 2015 34.747174 35.353766
## May 2015 46.800871 47.376541
## Jun 2015 39.557334 39.571802
## Jul 2015 34.053485 33.573768
## Aug 2015 36.228435 35.833488
## Sep 2015 24.918635 24.478713
## Oct 2015 8.672086 8.605752
## Nov 2015 -40.677195 -41.064456
## Dec 2015 -94.908579 -94.921414
## Jan 2016 -55.383135 -55.125978
## Feb 2016 -29.924007 -29.349144
## Mar 2016 -4.072846 -3.878576
## Apr 2016 34.687408 35.230024
## May 2016 46.785793 47.131999
## Jun 2016 39.513737 39.417012
## Jul 2016 34.200448 33.698854
## Aug 2016 36.169383 35.964420
## Sep 2016 NA NA
## Oct 2016 NA NA
## Nov 2016 NA NA
## Dec 2016 NA NA
## Jan 2017 NA NA
## Feb 2017 NA NA
##
## $figure
## [1] 3.998767 -34.759863 -46.821312 -39.493936 -33.856298 -36.175662
## [7] -24.854130 -8.723305 40.661334 94.845043 55.312609 29.866751
##
## $type
## [1] "additive"
##
## attr(,"class")
## [1] "decomposed.ts"plot(add_decomp$x[,1:6])
plot(add_decomp$x[,7:11])
plot(add_decomp$seasonal[1:60], type = "l")
plot(add_decomp$trend[,1:6])
plot(add_decomp$trend[,7:11])
plot(add_decomp$random[,1:6])
plot(add_decomp$random[,7:11])
# multiplicative decomposition
multi_decomp_PM25 <- decompose(ts_by_month[,'PM2.5'], type = 'multiplicative')
multi_decomp <- decompose(ts_by_month, type = 'multiplicative')
multi_decomp_PM25## $x
## Jan Feb Mar Apr May Jun Jul
## 2013 106.69892 78.82778 81.71640 105.67153 68.14651
## 2014 113.49194 156.06399 98.66667 86.64444 59.18804 58.03889 85.52151
## 2015 103.28199 102.95833 87.32258 71.85556 55.34543 59.44583 61.57124
## 2016 75.26210 48.35201 95.18952 67.05417 55.25806 64.08889 74.10753
## 2017 132.08333 77.61905
## Aug Sep Oct Nov Dec
## 2013 59.41532 61.74306 98.97849 84.30972 92.30914
## 2014 64.41331 70.99903 114.25806 103.31708 67.94247
## 2015 45.68548 53.15000 77.78226 125.42639 165.81989
## 2016 50.55242 59.25139 83.71640 106.37361 156.83468
## 2017
##
## $seasonal
## Jan Feb Mar Apr May Jun Jul
## 2013 1.1536745 0.9261272 0.7002633 0.7396729 0.8814532
## 2014 1.1888050 1.2376261 1.1536745 0.9261272 0.7002633 0.7396729 0.8814532
## 2015 1.1888050 1.2376261 1.1536745 0.9261272 0.7002633 0.7396729 0.8814532
## 2016 1.1888050 1.2376261 1.1536745 0.9261272 0.7002633 0.7396729 0.8814532
## 2017 1.1888050 1.2376261
## Aug Sep Oct Nov Dec
## 2013 0.6454082 0.7402740 1.1617698 1.2774791 1.3474468
## 2014 0.6454082 0.7402740 1.1617698 1.2774791 1.3474468
## 2015 0.6454082 0.7402740 1.1617698 1.2774791 1.3474468
## 2016 0.6454082 0.7402740 1.1617698 1.2774791 1.3474468
## 2017
##
## $trend
## Jan Feb Mar Apr May Jun Jul Aug
## 2013 NA NA NA NA NA NA
## 2014 87.14031 88.07252 88.66643 89.68875 91.11737 90.89406 89.45337 86.81522
## 2015 81.22384 79.44558 77.92155 75.65801 75.05941 80.05886 82.96959 79.52683
## 2016 78.40901 79.13415 79.59116 80.09264 79.54603 78.37778 80.37095 83.95796
## 2017 NA NA
## Sep Oct Nov Dec
## 2013 91.94639 91.93741 91.32442 88.40104
## 2014 84.12981 83.04094 82.26463 82.16314
## 2015 77.57935 77.70708 77.50339 77.69321
## 2016 NA NA NA NA
## 2017
##
## $random
## Jan Feb Mar Apr May Jun Jul
## 2013 NA NA NA NA NA
## 2014 1.0955578 1.4317686 0.9645571 1.0431149 0.9276228 0.8632645 1.0846244
## 2015 1.0696223 1.0471341 0.9713723 1.0254979 1.0529682 1.0038581 0.8418984
## 2016 0.8074204 0.4936978 1.0366711 0.9039877 0.9920094 1.1054779 1.0460777
## 2017 NA NA
## Aug Sep Oct Nov Dec
## 2013 NA 0.9071121 0.9266773 0.7226649 0.7749536
## 2014 1.1495959 1.1400135 1.1843349 0.9831170 0.6136951
## 2015 0.8900821 0.9254749 0.8615884 1.2668186 1.5839518
## 2016 0.9329225 NA NA NA NA
## 2017
##
## $figure
## [1] 1.1536745 0.9261272 0.7002633 0.7396729 0.8814532 0.6454082 0.7402740
## [8] 1.1617698 1.2774791 1.3474468 1.1888050 1.2376261
##
## $type
## [1] "multiplicative"
##
## attr(,"class")
## [1] "decomposed.ts"plot(multi_decomp_PM25)
multi_decomp## $x
## PM2.5 PM10 SO2 NO2 CO O3
## Mar 2013 106.69892 130.19355 41.782387 63.42723 1563.6962 49.67568
## Apr 2013 78.82778 93.94167 22.014413 42.76404 1000.7417 50.78922
## May 2013 81.71640 129.03360 26.635918 52.47490 1094.1250 71.12242
## Jun 2013 105.67153 118.43681 14.815358 52.75863 1628.0139 69.63165
## Jul 2013 68.14651 73.75403 7.801311 45.31321 1223.7702 84.50378
## Aug 2013 59.41532 81.54301 6.996331 45.75652 1075.8065 86.94474
## Sep 2013 61.74306 117.49028 13.512820 57.73660 1326.5194 49.56147
## Oct 2013 98.97849 133.95833 15.550846 73.50200 1392.0699 28.22625
## Nov 2013 84.30972 111.67361 23.508333 58.00556 1575.4167 19.07793
## Dec 2013 92.30914 121.84812 44.816955 66.48069 2213.8441 12.87716
## Jan 2014 113.49194 152.32863 61.379032 79.61559 2171.3710 12.22043
## Feb 2014 156.06399 176.36310 67.059524 74.99405 2458.0357 18.06548
## Mar 2014 98.66667 157.95161 40.067204 66.39785 1489.1129 47.58333
## Apr 2014 86.64444 141.04306 18.758333 54.85417 1075.8333 68.79861
## May 2014 59.18804 112.96250 14.621640 47.10403 934.5430 92.11142
## Jun 2014 58.03889 81.15000 6.341667 46.40833 990.5556 103.63750
## Jul 2014 85.52151 110.45161 5.244624 40.52151 907.9301 104.85215
## Aug 2014 64.41331 91.18616 4.044624 49.18696 1027.1505 94.55497
## Sep 2014 70.99903 97.71708 5.679167 59.13194 1162.5000 53.78361
## Oct 2014 114.25806 143.46505 8.204301 76.72177 1372.3118 26.74059
## Nov 2014 103.31708 142.29014 14.245139 77.05292 1785.4167 22.32083
## Dec 2014 67.94247 111.75591 33.370027 63.90403 1873.9247 29.50847
## Jan 2015 103.28199 128.48925 39.075672 76.75376 2011.6935 24.22755
## Feb 2015 102.95833 133.78869 29.120536 59.45833 1557.1429 38.70536
## Mar 2015 87.32258 155.15323 23.263441 59.77419 1257.3925 54.28763
## Apr 2015 71.85556 121.17083 9.693056 49.62222 835.8333 78.03472
## May 2015 55.34543 93.78898 9.059140 41.90726 779.1667 98.06855
## Jun 2015 59.44583 81.00139 6.877778 37.44861 865.8333 92.75417
## Jul 2015 61.57124 77.40188 4.935484 38.15995 978.3602 95.05780
## Aug 2015 45.68548 67.98118 3.108871 35.40726 909.6774 100.36694
## Sep 2015 53.15000 68.63472 4.965278 43.25139 956.6667 58.35278
## Oct 2015 77.78226 98.78629 5.424731 55.75672 1043.9516 36.92473
## Nov 2015 125.42639 131.26806 14.548611 60.19444 2085.4167 14.15139
## Dec 2015 165.81989 175.36425 24.176075 73.86290 3024.3280 17.96102
## Jan 2016 75.26210 88.44892 24.372312 52.94892 1639.1129 28.16935
## Feb 2016 48.35201 61.49856 17.149425 34.41810 933.9080 47.92385
## Mar 2016 95.18952 130.83737 20.548387 53.99194 1223.9247 48.20565
## Apr 2016 67.05417 109.42500 11.688889 41.57083 852.0833 72.54167
## May 2016 55.25806 95.65995 9.151882 37.45833 724.3280 99.22043
## Jun 2016 64.08889 82.17778 5.462500 34.99444 865.8333 126.60556
## Jul 2016 74.10753 82.06452 3.322581 34.48925 997.9839 96.86290
## Aug 2016 50.55242 64.67608 2.659946 39.72849 932.6613 76.13575
## Sep 2016 59.25139 70.65417 2.827778 50.75833 870.1389 56.39444
## Oct 2016 83.71640 94.86962 3.310484 61.01210 1246.5054 20.25538
## Nov 2016 106.37361 136.37778 10.468056 71.47778 1812.7778 14.05556
## Dec 2016 156.83468 169.18011 17.418011 89.91532 2604.1667 15.15188
## Jan 2017 132.08333 149.47984 18.594086 71.23522 2254.3011 27.93952
## Feb 2017 77.61905 90.79464 18.520833 56.49405 1173.6607 40.98810
## TEMP PRES DEWP RAIN WSPM
## Mar 2013 6.373252688 1011.6395 -6.3188172 0.0236559140 1.753495
## Apr 2013 12.599166667 1007.2869 -2.8651389 0.0141666667 2.223472
## May 2013 21.787903226 1002.1239 8.4520161 0.0076612903 1.771237
## Jun 2013 23.783055556 999.0772 17.5343056 0.1015277778 1.355000
## Jul 2013 27.156451613 995.0235 21.1831989 0.2802419355 1.359946
## Aug 2013 27.257661290 998.2694 20.0935484 0.0907258065 1.416398
## Sep 2013 20.487361111 1008.3138 14.5504167 0.1150000000 1.160972
## Oct 2013 13.199865591 1015.2367 5.3115591 0.0170698925 1.337366
## Nov 2013 5.880555556 1016.0646 -7.1959722 0.0008333333 1.769028
## Dec 2013 -0.005241935 1019.1534 -13.2763441 0.0000000000 1.699866
## Jan 2014 0.173655914 1019.2828 -12.6904570 0.0000000000 1.618548
## Feb 2014 0.089434524 1021.6240 -9.5802083 0.0098214286 1.438244
## Mar 2014 10.513844086 1012.8694 -5.8724462 0.0001344086 1.810215
## Apr 2014 17.184027778 1009.5782 3.9636111 0.0180555556 1.615556
## May 2014 21.865188172 1000.7228 6.9375000 0.0915322581 2.020027
## Jun 2014 24.901944444 999.3411 16.3420833 0.1980555556 1.393472
## Jul 2014 28.287634409 997.8254 19.5525538 0.0581989247 1.479973
## Aug 2014 26.208736559 1001.8782 18.1295699 0.1381720430 1.395833
## Sep 2014 21.123472222 1007.9142 13.8113889 0.0547222222 1.310833
## Oct 2014 14.044892473 1013.6585 6.2706989 0.0182795699 1.217070
## Nov 2014 6.371111111 1021.8025 -5.1972222 0.0004166667 1.538472
## Dec 2014 -0.226478495 1025.9004 -16.1875000 0.0000000000 2.189113
## Jan 2015 -0.125537634 1020.9488 -13.7850806 0.0005376344 1.604704
## Feb 2015 1.743154762 1017.6015 -12.4601190 0.0168154762 1.811905
## Mar 2015 8.439861751 1016.8734 -9.8809140 0.0088709677 2.213575
## Apr 2015 15.630138889 1010.8733 1.6945833 0.0556944444 2.483056
## May 2015 21.455510753 1003.4489 6.6134409 0.0501344086 2.471102
## Jun 2015 24.537777778 1000.3461 14.3013889 0.1362500000 1.980278
## Jul 2015 26.397849462 1001.2784 18.5057796 0.2701612903 1.627285
## Aug 2015 26.435752688 1003.4042 18.4461022 0.1258064516 1.615457
## Sep 2015 20.600833333 1011.3549 13.7698611 0.1372222222 1.601944
## Oct 2015 14.319758065 1015.7254 4.1530914 0.0190860215 1.848790
## Nov 2015 2.847777778 1023.4250 -1.2023611 0.0522222222 1.548056
## Dec 2015 -0.362768817 1024.0706 -7.2271505 0.0028225806 1.856452
## Jan 2016 -4.550000000 1025.9500 -18.0596774 0.0006720430 2.115054
## Feb 2016 1.398419540 1022.5598 -15.0882184 0.0181034483 2.449138
## Mar 2016 9.201881720 1016.7401 -8.7526882 0.0000000000 2.049866
## Apr 2016 16.454722222 1007.8117 -0.3329167 0.0066666667 2.382361
## May 2016 21.672849462 1005.5414 6.2018817 0.0458333333 2.207796
## Jun 2016 25.735833333 1000.7619 14.6834722 0.0991666667 1.803056
## Jul 2016 27.699462366 1000.6353 20.7272849 0.4233870968 1.698925
## Aug 2016 27.863306452 1003.5901 19.2868280 0.0724462366 1.618414
## Sep 2016 22.197344026 1009.3406 14.0559722 0.1598611111 1.554444
## Oct 2016 13.338648073 1017.1267 7.2780914 0.1174731183 1.476882
## Nov 2016 4.414027778 1020.5669 -3.7584722 0.0073611111 1.557778
## Dec 2016 0.640725806 1023.4661 -9.2724462 0.0000000000 1.435618
## Jan 2017 -1.140412186 1025.0380 -12.7756720 0.0004032258 1.969220
## Feb 2017 2.661532738 1022.0206 -13.6276786 0.0061011905 1.953869
##
## $seasonal
## Jan Feb Mar Apr May Jun Jul
## 2013 0.5789615 0.9512677 1.2078615 1.6021536 1.8483973
## 2014 0.2785447 0.3692298 0.5789615 0.9512677 1.2078615 1.6021536 1.8483973
## 2015 0.2785447 0.3692298 0.5789615 0.9512677 1.2078615 1.6021536 1.8483973
## 2016 0.2785447 0.3692298 0.5789615 0.9512677 1.2078615 1.6021536 1.8483973
## 2017 0.2785447 0.3692298 0.5789615 0.9512677 1.2078615 1.6021536 1.8483973
## 2018 0.2785447 0.3692298 0.5789615 0.9512677 1.2078615 1.6021536 1.8483973
## 2019 0.2785447 0.3692298 0.5789615 0.9512677 1.2078615 1.6021536 1.8483973
## 2020 0.2785447 0.3692298 0.5789615 0.9512677 1.2078615 1.6021536 1.8483973
## 2021 0.2785447 0.3692298 0.5789615 0.9512677 1.2078615 1.6021536 1.8483973
## 2022 0.2785447 0.3692298 0.5789615 0.9512677 1.2078615 1.6021536 1.8483973
## 2023 0.2785447 0.3692298 0.5789615 0.9512677 1.2078615 1.6021536 1.8483973
## 2024 0.2785447 0.3692298 0.5789615 0.9512677 1.2078615 1.6021536 1.8483973
## 2025 0.2785447 0.3692298 0.5789615 0.9512677 1.2078615 1.6021536 1.8483973
## 2026 0.2785447 0.3692298 0.5789615 0.9512677 1.2078615 1.6021536 1.8483973
## 2027 0.2785447 0.3692298 0.5789615 0.9512677 1.2078615 1.6021536 1.8483973
## 2028 0.2785447 0.3692298 0.5789615 0.9512677 1.2078615 1.6021536 1.8483973
## 2029 0.2785447 0.3692298 0.5789615 0.9512677 1.2078615 1.6021536 1.8483973
## 2030 0.2785447 0.3692298 0.5789615 0.9512677 1.2078615 1.6021536 1.8483973
## 2031 0.2785447 0.3692298 0.5789615 0.9512677 1.2078615 1.6021536 1.8483973
## 2032 0.2785447 0.3692298 0.5789615 0.9512677 1.2078615 1.6021536 1.8483973
## 2033 0.2785447 0.3692298 0.5789615 0.9512677 1.2078615 1.6021536 1.8483973
## 2034 0.2785447 0.3692298 0.5789615 0.9512677 1.2078615 1.6021536 1.8483973
## 2035 0.2785447 0.3692298 0.5789615 0.9512677 1.2078615 1.6021536 1.8483973
## 2036 0.2785447 0.3692298 0.5789615 0.9512677 1.2078615 1.6021536 1.8483973
## 2037 0.2785447 0.3692298 0.5789615 0.9512677 1.2078615 1.6021536 1.8483973
## 2038 0.2785447 0.3692298 0.5789615 0.9512677 1.2078615 1.6021536 1.8483973
## 2039 0.2785447 0.3692298 0.5789615 0.9512677 1.2078615 1.6021536 1.8483973
## 2040 0.2785447 0.3692298 0.5789615 0.9512677 1.2078615 1.6021536 1.8483973
## 2041 0.2785447 0.3692298 0.5789615 0.9512677 1.2078615 1.6021536 1.8483973
## 2042 0.2785447 0.3692298 0.5789615 0.9512677 1.2078615 1.6021536 1.8483973
## 2043 0.2785447 0.3692298 0.5789615 0.9512677 1.2078615 1.6021536 1.8483973
## 2044 0.2785447 0.3692298 0.5789615 0.9512677 1.2078615 1.6021536 1.8483973
## 2045 0.2785447 0.3692298 0.5789615 0.9512677 1.2078615 1.6021536 1.8483973
## 2046 0.2785447 0.3692298 0.5789615 0.9512677 1.2078615 1.6021536 1.8483973
## 2047 0.2785447 0.3692298 0.5789615 0.9512677 1.2078615 1.6021536 1.8483973
## 2048 0.2785447 0.3692298 0.5789615 0.9512677 1.2078615 1.6021536 1.8483973
## 2049 0.2785447 0.3692298 0.5789615 0.9512677 1.2078615 1.6021536 1.8483973
## 2050 0.2785447 0.3692298 0.5789615 0.9512677 1.2078615 1.6021536 1.8483973
## 2051 0.2785447 0.3692298 0.5789615 0.9512677 1.2078615 1.6021536 1.8483973
## 2052 0.2785447 0.3692298 0.5789615 0.9512677 1.2078615 1.6021536 1.8483973
## 2053 0.2785447 0.3692298 0.5789615 0.9512677 1.2078615 1.6021536 1.8483973
## 2054 0.2785447 0.3692298 0.5789615 0.9512677 1.2078615 1.6021536 1.8483973
## 2055 0.2785447 0.3692298 0.5789615 0.9512677 1.2078615 1.6021536 1.8483973
## 2056 0.2785447 0.3692298 0.5789615 0.9512677 1.2078615 1.6021536 1.8483973
## 2057 0.2785447 0.3692298
## Aug Sep Oct Nov Dec
## 2013 1.6237867 1.4283454 1.0022651 0.6628109 0.4463759
## 2014 1.6237867 1.4283454 1.0022651 0.6628109 0.4463759
## 2015 1.6237867 1.4283454 1.0022651 0.6628109 0.4463759
## 2016 1.6237867 1.4283454 1.0022651 0.6628109 0.4463759
## 2017 1.6237867 1.4283454 1.0022651 0.6628109 0.4463759
## 2018 1.6237867 1.4283454 1.0022651 0.6628109 0.4463759
## 2019 1.6237867 1.4283454 1.0022651 0.6628109 0.4463759
## 2020 1.6237867 1.4283454 1.0022651 0.6628109 0.4463759
## 2021 1.6237867 1.4283454 1.0022651 0.6628109 0.4463759
## 2022 1.6237867 1.4283454 1.0022651 0.6628109 0.4463759
## 2023 1.6237867 1.4283454 1.0022651 0.6628109 0.4463759
## 2024 1.6237867 1.4283454 1.0022651 0.6628109 0.4463759
## 2025 1.6237867 1.4283454 1.0022651 0.6628109 0.4463759
## 2026 1.6237867 1.4283454 1.0022651 0.6628109 0.4463759
## 2027 1.6237867 1.4283454 1.0022651 0.6628109 0.4463759
## 2028 1.6237867 1.4283454 1.0022651 0.6628109 0.4463759
## 2029 1.6237867 1.4283454 1.0022651 0.6628109 0.4463759
## 2030 1.6237867 1.4283454 1.0022651 0.6628109 0.4463759
## 2031 1.6237867 1.4283454 1.0022651 0.6628109 0.4463759
## 2032 1.6237867 1.4283454 1.0022651 0.6628109 0.4463759
## 2033 1.6237867 1.4283454 1.0022651 0.6628109 0.4463759
## 2034 1.6237867 1.4283454 1.0022651 0.6628109 0.4463759
## 2035 1.6237867 1.4283454 1.0022651 0.6628109 0.4463759
## 2036 1.6237867 1.4283454 1.0022651 0.6628109 0.4463759
## 2037 1.6237867 1.4283454 1.0022651 0.6628109 0.4463759
## 2038 1.6237867 1.4283454 1.0022651 0.6628109 0.4463759
## 2039 1.6237867 1.4283454 1.0022651 0.6628109 0.4463759
## 2040 1.6237867 1.4283454 1.0022651 0.6628109 0.4463759
## 2041 1.6237867 1.4283454 1.0022651 0.6628109 0.4463759
## 2042 1.6237867 1.4283454 1.0022651 0.6628109 0.4463759
## 2043 1.6237867 1.4283454 1.0022651 0.6628109 0.4463759
## 2044 1.6237867 1.4283454 1.0022651 0.6628109 0.4463759
## 2045 1.6237867 1.4283454 1.0022651 0.6628109 0.4463759
## 2046 1.6237867 1.4283454 1.0022651 0.6628109 0.4463759
## 2047 1.6237867 1.4283454 1.0022651 0.6628109 0.4463759
## 2048 1.6237867 1.4283454 1.0022651 0.6628109 0.4463759
## 2049 1.6237867 1.4283454 1.0022651 0.6628109 0.4463759
## 2050 1.6237867 1.4283454 1.0022651 0.6628109 0.4463759
## 2051 1.6237867 1.4283454 1.0022651 0.6628109 0.4463759
## 2052 1.6237867 1.4283454 1.0022651 0.6628109 0.4463759
## 2053 1.6237867 1.4283454 1.0022651 0.6628109 0.4463759
## 2054 1.6237867 1.4283454 1.0022651 0.6628109 0.4463759
## 2055 1.6237867 1.4283454 1.0022651 0.6628109 0.4463759
## 2056 1.6237867 1.4283454 1.0022651 0.6628109 0.4463759
## 2057
##
## $trend
## [,1] [,2] [,3] [,4] [,5] [,6] [,7]
## Mar 2013 NA NA NA NA NA NA NA
## Apr 2013 NA NA NA NA NA NA NA
## May 2013 NA NA NA NA NA NA NA
## Jun 2013 NA NA NA NA NA NA NA
## Jul 2013 NA NA NA NA NA NA NA
## Aug 2013 NA NA NA NA NA NA NA
## Sep 2013 91.94639 121.20365 28.75130 59.52619 1557.177 45.97084 13.40445
## Oct 2013 91.93741 124.32279 28.54417 60.15372 1557.198 46.63405 13.76801
## Nov 2013 91.32442 125.61572 27.90790 60.43369 1553.677 48.25898 13.96227
## Dec 2013 88.40104 123.39247 27.05424 59.94531 1520.467 50.55043 14.01211
## Jan 2014 87.14031 123.36792 26.59464 59.48106 1480.747 52.81519 14.10586
## Feb 2014 88.07252 125.29879 26.36512 59.42434 1465.559 53.98013 14.10929
## Mar 2014 88.66643 124.87670 25.91573 59.62542 1456.698 54.47315 14.09209
## Apr 2014 89.68875 124.44893 25.28322 59.81771 1449.040 54.58717 14.15380
## May 2014 91.11737 126.12073 24.59115 60.74551 1456.967 54.66039 14.20945
## Jun 2014 90.89406 126.97591 23.72823 61.43179 1451.554 55.48848 14.22067
## Jul 2014 89.45337 125.56210 22.32197 61.20519 1430.737 56.68175 14.19899
## Aug 2014 86.81522 122.79486 19.81187 60.43862 1386.547 58.04204 14.25543
## Sep 2014 84.12981 120.90432 17.53092 59.51532 1339.355 59.18138 14.23792
## Oct 2014 83.04094 119.95971 16.45304 59.02133 1319.700 59.84556 14.08676
## Nov 2014 82.26463 118.33281 15.84355 58.58680 1303.226 60.47862 14.00494
## Dec 2014 82.16314 117.52772 15.63412 57.99695 1291.555 60.27336 13.97270
## Jan 2015 81.22384 116.14445 15.64358 57.52523 1289.293 59.41179 13.87878
## Feb 2015 79.44558 113.80051 15.59171 56.85268 1287.332 59.24585 13.80950
## Mar 2015 77.92155 111.62187 15.52297 55.61683 1273.861 59.67840 13.79718
## Apr 2015 75.65801 108.54849 15.37741 54.08160 1251.603 60.29312 13.78686
## May 2015 75.05941 106.22762 15.27424 52.50562 1250.422 60.37707 13.65151
## Jun 2015 80.05886 108.41871 14.90380 52.21813 1310.855 59.55553 13.49902
## Jul 2015 82.96959 109.40072 13.90808 51.64122 1343.264 59.23863 13.30899
## Aug 2015 79.52683 104.72028 12.79665 49.60601 1301.772 59.78697 13.11027
## Sep 2015 77.57935 100.69503 12.18472 48.32174 1274.409 59.91766 13.12766
## Oct 2015 77.70708 99.19246 12.15476 47.74534 1273.692 59.43537 13.19377
## Nov 2015 77.50339 98.78101 12.24178 47.22449 1272.084 59.25449 13.23718
## Dec 2015 77.69321 98.90798 12.18667 46.93686 1269.799 60.71295 13.29616
## Jan 2016 78.40901 99.15127 12.06050 46.68166 1270.617 62.19864 13.40031
## Feb 2016 79.13415 99.20784 11.97459 46.70876 1272.392 61.26422 13.51402
## Mar 2016 79.59116 99.15427 11.86682 47.20160 1269.745 60.17299 13.64003
## Apr 2016 80.09264 99.07522 11.68967 47.73337 1274.579 59.39684 13.66567
## May 2016 79.54603 99.12493 11.43155 48.42248 1271.659 58.69829 13.69005
## Jun 2016 78.37778 99.08016 10.97994 49.56147 1242.792 58.57725 13.79712
## Jul 2016 80.37095 101.36544 10.45759 50.99225 1250.918 58.45062 13.98100
## Aug 2016 83.95796 105.12907 10.27398 52.67401 1286.541 58.15206 14.17570
## Sep 2016 NA NA NA NA NA NA NA
## Oct 2016 NA NA NA NA NA NA NA
## Nov 2016 NA NA NA NA NA NA NA
## Dec 2016 NA NA NA NA NA NA NA
## Jan 2017 NA NA NA NA NA NA NA
## Feb 2017 NA NA NA NA NA NA NA
## [,8] [,9] [,10] [,11]
## Mar 2013 NA NA NA NA
## Apr 2013 NA NA NA NA
## May 2013 NA NA NA NA
## Jun 2013 NA NA NA NA
## Jul 2013 NA NA NA NA
## Aug 2013 NA NA NA NA
## Sep 2013 1009.476 2.951774 0.05407861 1.577661
## Oct 2013 1009.623 3.254904 0.05326058 1.554694
## Nov 2013 1009.660 3.476331 0.05691724 1.539731
## Dec 2013 1009.612 3.363550 0.06443386 1.551700
## Jan 2014 1009.740 3.245931 0.05920406 1.558304
## Feb 2014 1010.007 3.096155 0.05192919 1.562449
## Mar 2014 1010.141 2.983529 0.05139454 1.567836
## Apr 2014 1010.058 2.992701 0.04893337 1.569068
## May 2014 1010.232 3.115946 0.04896641 1.554449
## Jun 2014 1010.752 3.077929 0.04894905 1.565228
## Jul 2014 1011.103 2.911022 0.04897145 1.585036
## Aug 2014 1011.004 2.745416 0.04928527 1.600029
## Sep 2014 1011.004 2.458400 0.04994072 1.632404
## Oct 2014 1011.224 2.196838 0.05187303 1.685357
## Nov 2014 1011.392 2.088793 0.05171640 1.740298
## Dec 2014 1011.547 1.990261 0.04741626 1.783543
## Jan 2015 1011.733 1.861617 0.05367280 1.814131
## Feb 2015 1011.941 1.831190 0.06198933 1.829420
## Mar 2015 1012.148 1.842649 0.06491159 1.850700
## Apr 2015 1012.377 1.752685 0.06838270 1.889152
## May 2015 1012.531 1.830903 0.07057486 1.915873
## Jun 2015 1012.522 2.370704 0.07285104 1.902411
## Jul 2015 1012.654 2.565944 0.07297424 1.909815
## Aug 2015 1013.069 2.278331 0.07303351 1.957631
## Sep 2015 1013.270 2.215837 0.07271755 1.977361
## Oct 2015 1013.137 2.178367 0.07030510 1.966344
## Nov 2015 1013.097 2.076739 0.06808307 1.951177
## Dec 2015 1013.201 2.075511 0.06635872 1.932822
## Jan 2016 1013.192 2.183994 0.07119799 1.928423
## Feb 2016 1013.173 2.311587 0.07535905 1.931531
## Mar 2016 1013.097 2.358539 0.07407900 1.929675
## Apr 2016 1013.071 2.500668 0.07912175 1.912200
## May 2016 1013.010 2.524372 0.08135200 1.897109
## Jun 2016 1012.866 2.332647 0.07936518 1.879979
## Jul 2016 1012.803 2.467593 0.07923637 1.856368
## Aug 2016 1012.742 2.748616 0.07872507 1.829655
## Sep 2016 NA NA NA NA
## Oct 2016 NA NA NA NA
## Nov 2016 NA NA NA NA
## Dec 2016 NA NA NA NA
## Jan 2017 NA NA NA NA
## Feb 2017 NA NA NA NA
##
## $random
## x/seasonal.x.PM2.5 x/seasonal.x.PM10 x/seasonal.x.SO2 x/seasonal.x.NO2
## Mar 2013 NA NA NA NA
## Apr 2013 NA NA NA NA
## May 2013 NA NA NA NA
## Jun 2013 NA NA NA NA
## Jul 2013 NA NA NA NA
## Aug 2013 NA NA NA NA
## Sep 2013 0.4701324 0.6786612 0.3290449 0.6790627
## Oct 2013 1.0741526 1.0750691 0.5435682 1.2191412
## Nov 2013 1.3928397 1.3412723 1.2708814 1.4481077
## Dec 2013 2.3393035 2.2122256 3.7111312 2.4845030
## Jan 2014 4.6757474 4.4328640 8.2857359 4.8053447
## Feb 2014 4.7991632 3.8120982 6.8886468 3.4179495
## Mar 2014 1.9220363 2.1847060 2.6703976 1.9234147
## Apr 2014 1.0155470 1.1914006 0.7799362 0.9640001
## May 2014 0.5377936 0.7415333 0.4922663 0.6419878
## Jun 2014 0.3985469 0.3988991 0.1668145 0.4715184
## Jul 2014 0.5172295 0.4759027 0.1271120 0.3581806
## Aug 2014 0.4569311 0.4573195 0.1257256 0.5011947
## Sep 2014 0.5908391 0.5658423 0.2268019 0.6956010
## Oct 2014 1.3728149 1.1932408 0.4975225 1.2969613
## Nov 2014 1.8948262 1.8141785 1.3565146 1.9842650
## Dec 2014 1.8525228 2.1302445 4.7817002 2.4684388
## Jan 2015 4.5650571 3.9716725 8.9675850 4.7901207
## Feb 2015 3.5099018 3.1840393 5.0583505 2.8324682
## Mar 2015 1.9356166 2.4008325 2.5885075 1.8563413
## Apr 2015 0.9983957 1.1734689 0.6626356 0.9645483
## May 2015 0.6104632 0.7309661 0.4910325 0.6607944
## Jun 2015 0.4634553 0.4663200 0.2880361 0.4476207
## Jul 2015 0.4014797 0.3827684 0.1919849 0.3997753
## Aug 2015 0.3537819 0.3997872 0.1496158 0.4395710
## Sep 2015 0.4796494 0.4772024 0.2852953 0.6266489
## Oct 2015 0.9987052 0.9936545 0.4452966 1.1651549
## Nov 2015 2.4416231 2.0049151 1.7930293 1.9230897
## Dec 2015 4.7813753 3.9719975 4.4442642 3.5254259
## Jan 2016 3.4460015 3.2025756 7.2549857 4.0720772
## Feb 2016 1.6548318 1.6788899 3.8787532 1.9956841
## Mar 2016 2.0657351 2.2791386 2.9908433 1.9757068
## Apr 2016 0.8800968 1.1610443 1.0511590 0.9155118
## May 2016 0.5751221 0.7989694 0.6628087 0.6404486
## Jun 2016 0.5103706 0.5176826 0.3105184 0.4407078
## Jul 2016 0.4988476 0.4379960 0.1718891 0.3659184
## Aug 2016 0.3708097 0.3788714 0.1594429 0.4644904
## Sep 2016 NA NA NA NA
## Oct 2016 NA NA NA NA
## Nov 2016 NA NA NA NA
## Dec 2016 NA NA NA NA
## Jan 2017 NA NA NA NA
## Feb 2017 NA NA NA NA
## x/seasonal.x.CO x/seasonal.x.O3 x/seasonal.x.TEMP x/seasonal.x.PRES
## Mar 2013 NA NA NA NA
## Apr 2013 NA NA NA NA
## May 2013 NA NA NA NA
## Jun 2013 NA NA NA NA
## Jul 2013 NA NA NA NA
## Aug 2013 NA NA NA NA
## Sep 2013 0.5964067 0.7547941 1.0700491783 0.6993048
## Oct 2013 0.8919380 0.6039034 0.9565675996 1.0032880
## Nov 2013 1.5298365 0.5964356 0.6354374035 1.5182968
## Dec 2013 3.2618892 0.5706826 -0.0008380838 2.2614352
## Jan 2014 5.2645158 0.8306781 0.0441972336 3.6240171
## Feb 2014 4.5424280 0.9063975 0.0171673513 2.7394908
## Mar 2014 1.7656659 1.5087689 1.2886545842 1.7318961
## Apr 2014 0.7804800 1.3249100 1.2762892222 1.0507290
## May 2014 0.5310464 1.3951589 1.2739687396 0.8201167
## Jun 2014 0.4259333 1.1657621 1.0929717397 0.6171134
## Jul 2014 0.3433185 1.0007804 1.0778141675 0.5339051
## Aug 2014 0.4562161 1.0032582 1.1322357224 0.6102853
## Sep 2014 0.6076649 0.6362556 1.0386892246 0.6979714
## Oct 2014 1.0375168 0.4458168 0.9947749536 1.0001416
## Nov 2014 2.0669520 0.5568253 0.6863479112 1.5242558
## Dec 2014 3.2504131 1.0967825 -0.0363116507 2.2720517
## Jan 2015 5.6016435 1.4640034 -0.0324733968 3.6227896
## Feb 2015 3.2759779 1.7693606 0.3418701821 2.7234910
## Mar 2015 1.7049005 1.5712095 1.0565627263 1.7352954
## Apr 2015 0.7020212 1.3605590 1.1917764352 1.0496674
## May 2015 0.5158896 1.3447471 1.3011914054 0.8204836
## Jun 2015 0.4122640 0.9720916 1.1345634475 0.6166541
## Jul 2015 0.3940415 0.8681353 1.0730701936 0.5349317
## Aug 2015 0.4303517 1.0338442 1.2417979821 0.6099690
## Sep 2015 0.5255553 0.6818258 1.0986623118 0.6987873
## Oct 2015 0.8177741 0.6198545 1.0828896712 1.0002889
## Nov 2015 2.4733603 0.3603199 0.3245793951 1.5241071
## Dec 2015 5.3357205 0.6627487 -0.0611227761 2.2642971
## Jan 2016 4.6312622 1.6259273 -1.2189943480 3.6352950
## Feb 2016 1.9878625 2.1185953 0.2802566940 2.7334329
## Mar 2016 1.6649021 1.3837150 1.1652302497 1.7334425
## Apr 2016 0.7027690 1.2838712 1.2657761794 1.0457714
## May 2016 0.4715715 1.3994537 1.3106714925 0.8218053
## Jun 2016 0.4348422 1.3490239 1.1642481136 0.6167009
## Jul 2016 0.4316177 0.8965469 1.0718592053 0.5345097
## Aug 2016 0.4464486 0.8062961 1.2104845181 0.6102789
## Sep 2016 NA NA NA NA
## Oct 2016 NA NA NA NA
## Nov 2016 NA NA NA NA
## Dec 2016 NA NA NA NA
## Jan 2017 NA NA NA NA
## Feb 2017 NA NA NA NA
## x/seasonal.x.DEWP x/seasonal.x.RAIN x/seasonal.x.WSPM
## Mar 2013 NA NA NA
## Apr 2013 NA NA NA
## May 2013 NA NA NA
## Jun 2013 NA NA NA
## Jul 2013 NA NA NA
## Aug 2013 NA NA NA
## Sep 2013 3.4511118 1.488809394 0.5151989
## Oct 2013 1.6281753 0.319773353 0.8582672
## Nov 2013 -3.1230488 0.022089470 1.7334058
## Dec 2013 -8.8425974 0.000000000 2.4541779
## Jan 2014 -14.0359936 0.000000000 3.7288813
## Feb 2014 -8.3802223 0.512231566 2.4930449
## Mar 2014 -3.3996882 0.004517107 1.9942515
## Apr 2014 1.3922750 0.387884998 1.0823743
## May 2014 1.8432996 1.547600117 1.0758793
## Jun 2014 3.3139400 2.525448835 0.5556696
## Jul 2014 3.6338142 0.642949165 0.5051488
## Aug 2014 4.0667768 1.726529543 0.5372505
## Sep 2014 3.9332496 0.767141930 0.5621943
## Oct 2014 2.8479693 0.351594261 0.7205116
## Nov 2014 -3.7539314 0.012155444 1.3337564
## Dec 2014 -18.2208623 0.000000000 2.7496911
## Jan 2015 -26.5842339 0.035961514 3.1756419
## Feb 2015 -18.4285853 0.734675356 2.6824103
## Mar 2015 -9.2620038 0.236047321 2.0658963
## Apr 2015 1.0163807 0.856175750 1.3817098
## May 2015 2.9905075 0.588123761 1.0678418
## Jun 2015 3.7652753 1.167337983 0.6497070
## Jul 2015 3.9017995 2.002895060 0.4609746
## Aug 2015 4.9860743 1.060844563 0.5082011
## Sep 2015 4.3506946 1.321149496 0.5671896
## Oct 2015 1.9022074 0.270860670 0.9380922
## Nov 2015 -0.8735008 1.157248481 1.1970166
## Dec 2015 -7.8008382 0.095290065 2.1517460
## Jan 2016 -29.6868191 0.033887106 3.9375337
## Feb 2016 -17.6779109 0.650622523 3.4341147
## Mar 2016 -6.4098642 0.000000000 1.8348120
## Apr 2016 -0.1399512 0.088574791 1.3096994
## May 2016 2.0340097 0.466440336 0.9634952
## Jun 2016 3.9289423 0.779886790 0.5986210
## Jul 2016 4.5443690 2.890797901 0.4951249
## Aug 2016 4.3213338 0.566726820 0.5447426
## Sep 2016 NA NA NA
## Oct 2016 NA NA NA
## Nov 2016 NA NA NA
## Dec 2016 NA NA NA
## Jan 2017 NA NA NA
## Feb 2017 NA NA NA
##
## $figure
## [1] 0.5789615 0.9512677 1.2078615 1.6021536 1.8483973 1.6237867 1.4283454
## [8] 1.0022651 0.6628109 0.4463759 0.2785447 0.3692298
##
## $type
## [1] "multiplicative"
##
## attr(,"class")
## [1] "decomposed.ts"plot(multi_decomp$x[,1:6])
plot(multi_decomp$x[,7:11])
plot(multi_decomp$seasonal[1:60], type = "l")
plot(multi_decomp$trend[,1:6])
plot(multi_decomp$trend[,7:11])
plot(multi_decomp$random[,1:6])
plot(multi_decomp$random[,7:11])
Moving average smoothing (monthly)
STL.PM2.5<- stl( ts_by_month[,'PM2.5'] , t.window = 13, s.window = 'periodic')
plot(STL.PM2.5)
STL.PM10<- stl( ts_by_month[,'PM10'] , t.window = 13, s.window = 'periodic')
plot(STL.PM10)
STL.NO2<- stl( ts_by_month[,'NO2'] , t.window = 13, s.window = 'periodic')
plot(STL.NO2)
STL.SO2<- stl( ts_by_month[,'NO2'] , t.window = 13, s.window = 'periodic')
plot(STL.SO2)
STL.CO<- stl( ts_by_month[,'CO'] , t.window = 13, s.window = 'periodic')
plot(STL.CO)
STL.O3<- stl( ts_by_month[,'O3'] , t.window = 13, s.window = 'periodic')
plot(STL.O3)
STL.TEMP<- stl( ts_by_month[,'TEMP'] , t.window = 13, s.window = 'periodic')
plot(STL.TEMP)
STL.PRES<- stl( ts_by_month[,'PRES'] , t.window = 13, s.window = 'periodic')
plot(STL.PRES)
STL.RAIN<- stl( ts_by_month[,'RAIN'] , t.window = 13, s.window = 'periodic')
plot(STL.RAIN)
Tesk for random / remainder using Ljung-Box Test (monthly)
auto.arima(remainder(STL.PM2.5))## Series: remainder(STL.PM2.5)
## ARIMA(0,0,0) with zero mean
##
## sigma^2 = 357.8: log likelihood = -209.23
## AIC=420.46 AICc=420.55 BIC=422.33x <- lapply(1:25, function(i){
p <- Box.test(remainder(STL.PM2.5), lag = i, type = "Ljung-Box")
output <- data.frame(lag = i, p_value = p$p.value)
return(output) }) %>% bind_rows
plot(x = x$lag,
y = x$p_value, ylim = c(0,1),
main = "Series white_noise (PM2.5) - Ljung-Box Test",
xlab = "Lag", ylab = "P-Value")
abline(h = 0.05, col="red", lwd=3, lty=2)
auto.arima(remainder(STL.PM10))## Series: remainder(STL.PM10)
## ARIMA(2,0,1) with zero mean
##
## Coefficients:
## ar1 ar2 ma1
## 0.5851 -0.6317 -0.8254
## s.e. 0.1170 0.1137 0.0924
##
## sigma^2 = 180: log likelihood = -192.29
## AIC=392.57 AICc=393.5 BIC=400.06x <- lapply(1:25, function(i){
p <- Box.test(remainder(STL.PM10), lag = i, type = "Ljung-Box")
output <- data.frame(lag = i, p_value = p$p.value)
return(output) }) %>% bind_rows
plot(x = x$lag,
y = x$p_value, ylim = c(0,1),
main = "Series white_noise(PM10) - Ljung-Box Test",
xlab = "Lag", ylab = "P-Value")
abline(h = 0.05, col="red", lwd=3, lty=2)
auto.arima(remainder(STL.SO2))## Series: remainder(STL.SO2)
## ARIMA(0,0,0) with zero mean
##
## sigma^2 = 32.73: log likelihood = -151.83
## AIC=305.65 AICc=305.74 BIC=307.52x <- lapply(1:25, function(i){
p <- Box.test(remainder(STL.SO2), lag = i, type = "Ljung-Box")
output <- data.frame(lag = i, p_value = p$p.value)
return(output) }) %>% bind_rows
plot(x = x$lag,
y = x$p_value, ylim = c(0,1),
main = "Series white_noise(SO2) - Ljung-Box Test",
xlab = "Lag", ylab = "P-Value")
abline(h = 0.05, col="red", lwd=3, lty=2)
auto.arima(remainder(STL.NO2))## Series: remainder(STL.NO2)
## ARIMA(0,0,0) with zero mean
##
## sigma^2 = 32.73: log likelihood = -151.83
## AIC=305.65 AICc=305.74 BIC=307.52x <- lapply(1:25, function(i){
p <- Box.test(remainder(STL.NO2), lag = i, type = "Ljung-Box")
output <- data.frame(lag = i, p_value = p$p.value)
return(output) }) %>% bind_rows
plot(x = x$lag,
y = x$p_value, ylim = c(0,1),
main = "Series white_noise(NO2) - Ljung-Box Test",
xlab = "Lag", ylab = "P-Value")
abline(h = 0.05, col="red", lwd=3, lty=2)
auto.arima(remainder(STL.CO))## Series: remainder(STL.CO)
## ARIMA(0,0,0) with zero mean
##
## sigma^2 = 53698: log likelihood = -329.5
## AIC=660.99 AICc=661.08 BIC=662.86x <- lapply(1:25, function(i){
p <- Box.test(remainder(STL.CO), lag = i, type = "Ljung-Box")
output <- data.frame(lag = i, p_value = p$p.value)
return(output) }) %>% bind_rows
plot(x = x$lag,
y = x$p_value, ylim = c(0,1),
main = "Series white_noise(CO) - Ljung-Box Test",
xlab = "Lag", ylab = "P-Value")
abline(h = 0.05, col="red", lwd=3, lty=2)
auto.arima(remainder(STL.O3))## Series: remainder(STL.O3)
## ARIMA(0,0,0) with zero mean
##
## sigma^2 = 51.6: log likelihood = -162.76
## AIC=327.51 AICc=327.6 BIC=329.38x <- lapply(1:25, function(i){
p <- Box.test(remainder(STL.O3), lag = i, type = "Ljung-Box")
output <- data.frame(lag = i, p_value = p$p.value)
return(output) }) %>% bind_rows
plot(x = x$lag,
y = x$p_value, ylim = c(0,1),
main = "Series white_noise(O3) - Ljung-Box Test",
xlab = "Lag", ylab = "P-Value")
abline(h = 0.05, col="red", lwd=3, lty=2)
auto.arima(remainder(STL.TEMP))## Series: remainder(STL.TEMP)
## ARIMA(1,0,1) with zero mean
##
## Coefficients:
## ar1 ma1
## 0.3815 -0.8530
## s.e. 0.2102 0.1376
##
## sigma^2 = 0.5325: log likelihood = -52.3
## AIC=110.6 AICc=111.14 BIC=116.21x <- lapply(1:25, function(i){
p <- Box.test(remainder(STL.TEMP), lag = i, type = "Ljung-Box")
output <- data.frame(lag = i, p_value = p$p.value)
return(output) }) %>% bind_rows
plot(x = x$lag,
y = x$p_value, ylim = c(0,1),
main = "Series white_noise(TEMP) - Ljung-Box Test",
xlab = "Lag", ylab = "P-Value")
abline(h = 0.05, col="red", lwd=3, lty=2)
auto.arima(remainder(STL.PRES))## Series: remainder(STL.PRES)
## ARIMA(0,0,0) with zero mean
##
## sigma^2 = 1.942: log likelihood = -84.04
## AIC=170.09 AICc=170.17 BIC=171.96x <- lapply(1:25, function(i){
p <- Box.test(remainder(STL.PRES), lag = i, type = "Ljung-Box")
output <- data.frame(lag = i, p_value = p$p.value)
return(output) }) %>% bind_rows
plot(x = x$lag,
y = x$p_value, ylim = c(0,1),
main = "Series white_noise(PRES) - Ljung-Box Test",
xlab = "Lag", ylab = "P-Value")
abline(h = 0.05, col="red", lwd=3, lty=2)
auto.arima(remainder(STL.RAIN))## Series: remainder(STL.RAIN)
## ARIMA(2,0,0) with zero mean
##
## Coefficients:
## ar1 ar2
## -0.5866 -0.2711
## s.e. 0.1379 0.1360
##
## sigma^2 = 0.001314: log likelihood = 91.96
## AIC=-177.91 AICc=-177.37 BIC=-172.3x <- lapply(1:25, function(i){
p <- Box.test(remainder(STL.RAIN), lag = i, type = "Ljung-Box")
output <- data.frame(lag = i, p_value = p$p.value)
return(output) }) %>% bind_rows
plot(x = x$lag,
y = x$p_value, ylim = c(0,1),
main = "Series white_noise(RAIN) - Ljung-Box Test",
xlab = "Lag", ylab = "P-Value")
abline(h = 0.05, col="red", lwd=3, lty=2)
forcesting for diffrent varible
forecast(STL.PM2.5, h=5, method = 'naive')## Point Forecast Lo 80 Hi 80 Lo 95 Hi 95
## Mar 2017 77.18490 44.111134 110.2587 26.60294 127.7669
## Apr 2017 56.62231 9.848936 103.3957 -14.91139 128.1560
## May 2017 43.71513 -13.570317 101.0006 -43.89540 131.3257
## Jun 2017 53.04516 -13.102381 119.1927 -48.11877 154.2091
## Jul 2017 53.96631 -19.988884 127.9215 -59.13840 167.0710STL.PM2.5 %>% forecast(method="naive", h=5) %>%
autoplot(ylab = 'PM2.5') + ggtitle('Forecasting')
forecast(STL.PM10, h=5, method = 'naive')## Point Forecast Lo 80 Hi 80 Lo 95 Hi 95
## Mar 2017 120.79043 88.517348 153.0635 71.43301 170.1479
## Apr 2017 93.15092 47.509882 138.7919 23.34898 162.9529
## May 2017 84.11630 28.217673 140.0149 -1.37327 169.6059
## Jun 2017 67.14942 2.603248 131.6956 -31.56543 165.8643
## Jul 2017 62.57884 -9.585969 134.7437 -47.78771 172.9454STL.PM10 %>% forecast(method="naive", h=5) %>%
autoplot(ylab = 'PM10') + ggtitle('Forecasting')
forecast(STL.SO2, h=5, method = 'naive')## Point Forecast Lo 80 Hi 80 Lo 95 Hi 95
## Mar 2017 61.73933 51.45773 72.02093 46.014975 77.46368
## Apr 2017 48.03199 33.49161 62.57237 25.794394 70.26958
## May 2017 45.55294 27.74468 63.36120 18.317562 72.78832
## Jun 2017 43.79880 23.23559 64.36200 12.350091 75.24750
## Jul 2017 40.59675 17.60639 63.58712 5.436031 75.75748STL.SO2 %>% forecast(method="naive", h=5) %>%
autoplot(ylab = 'SO2') + ggtitle('Forecasting')
forecast(STL.NO2, h=5, method = 'naive')## Point Forecast Lo 80 Hi 80 Lo 95 Hi 95
## Mar 2017 61.73933 51.45773 72.02093 46.014975 77.46368
## Apr 2017 48.03199 33.49161 62.57237 25.794394 70.26958
## May 2017 45.55294 27.74468 63.36120 18.317562 72.78832
## Jun 2017 43.79880 23.23559 64.36200 12.350091 75.24750
## Jul 2017 40.59675 17.60639 63.58712 5.436031 75.75748STL.NO2 %>% forecast(method="naive", h=5) %>%
autoplot(ylab = 'NO2') + ggtitle('Forecasting')
forecast(STL.CO, h=5, method = 'naive')## Point Forecast Lo 80 Hi 80 Lo 95 Hi 95
## Mar 2017 973.1480 572.96264 1373.333 361.1173 1585.179
## Apr 2017 541.0066 -24.94091 1106.954 -324.5354 1406.549
## May 2017 493.1916 -199.94972 1186.333 -566.8766 1553.260
## Jun 2017 705.9207 -94.44995 1506.291 -518.1406 1929.982
## Jul 2017 653.5837 -241.25789 1548.425 -714.9584 2022.126STL.CO %>% forecast(method="naive", h=5) %>%
autoplot(ylab = 'CO') + ggtitle('Forecasting')
forecast(STL.O3, h=5, method = 'naive')## Point Forecast Lo 80 Hi 80 Lo 95 Hi 95
## Mar 2017 57.51380 45.54433 69.48328 39.20807 75.81954
## Apr 2017 74.45243 57.52503 91.37983 48.56420 100.34065
## May 2017 96.37771 75.64596 117.10945 64.67124 128.08418
## Jun 2017 104.09568 80.15672 128.03464 67.48420 140.70716
## Jul 2017 100.94909 74.18452 127.71366 60.01621 141.88196STL.O3 %>% forecast(method="naive", h=5) %>%
autoplot(ylab = 'O3') + ggtitle('Forecasting')
forecast(STL.TEMP, h=5, method = 'naive')## Point Forecast Lo 80 Hi 80 Lo 95 Hi 95
## Mar 2017 10.39841 8.668394 12.12842 7.752579 13.04424
## Apr 2017 17.13377 14.687157 19.58038 13.391999 20.87554
## May 2017 23.26267 20.266196 26.25914 18.679959 27.84538
## Jun 2017 26.24717 22.787142 29.70720 20.955513 31.53883
## Jul 2017 28.83308 24.964650 32.70151 22.916826 34.74934STL.TEMP %>% forecast(method="naive", h=5) %>%
autoplot(ylab = 'TEMP') + ggtitle('Forecasting')
forecast(STL.PRES, h=5, method = 'naive')## Point Forecast Lo 80 Hi 80 Lo 95 Hi 95
## Mar 2017 1016.220 1013.4819 1018.959 1012.0323 1020.408
## Apr 2017 1010.481 1006.6084 1014.354 1004.5584 1016.404
## May 2017 1004.456 999.7135 1009.199 997.2028 1011.710
## Jun 2017 1001.340 995.8635 1006.817 992.9643 1009.716
## Jul 2017 1000.111 993.9874 1006.234 990.7461 1009.475STL.PRES %>% forecast(method="naive", h=5) %>%
autoplot(ylab = 'PRES') + ggtitle('Forecasting')
forecast(STL.RAIN, h=5, method = 'naive')## Point Forecast Lo 80 Hi 80 Lo 95 Hi 95
## Mar 2017 0.001549438 -0.09152873 0.09462761 -0.14080134 0.1439002
## Apr 2017 0.016914377 -0.11471803 0.14854679 -0.18440002 0.2182288
## May 2017 0.041943287 -0.11927283 0.20315941 -0.20461549 0.2885021
## Jun 2017 0.126730572 -0.05942577 0.31288691 -0.15797098 0.4114321
## Jul 2017 0.250805358 0.04267624 0.45893448 -0.06750066 0.5691114STL.RAIN %>% forecast(method="naive", h=5) %>%
autoplot(ylab = 'RAIN') + ggtitle('Forecasting')
Finding:
The air polution condition are correlated with thr weather condition.
Diffrence season has